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ABSTRACT 



The spectral analysis technique is used in order to predict the 
wave induced rolling motions of a ship in a seaway. Non-linear 
equations of motion for the ship, the anti-rolling tank, and 
for wave excitations are developed. Solution techniques for 
the equations are examined, and a linearized form is adopted. 

A design method using a digital computer program based upon 
numerical solution of the linearized equations is proposed and 
examined for a number of different representative ships. The 
design method is demonstrated to be capable for use to design 
passive anti-rolling tanks for ships and to determine feasi- 
bility of activated anti-rolling tanks. 
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I. Introduction 
1. General Discussion 

Successful roll stabilization of ships by means of 
fluid filled Anti-Rolling tanks requires a thorough analysis 
of the forces and moments involved. The author undertook 
this task with a goal of being able to use the results as a 
design tool of the naval architect. As such, the Anti-Rolling 
tank design method derived in this paper is intended to be 
a single step in the larger process of designing the whole 
ship. 

Modeling the dynamics of a ship with an Anti-Rolling 
tank has been investigated by many authors of whom only a 
small percentage are represented in the reference list. Sev- 
eral problems have prevented their analyses from being complete 
or from possessing sufficient accuracy in the range of interest. 

One central problem is that rolling is a resonance 
phenomenon. The ship is an underdamped dynamic system in roll, 
one which magnifies wave excitations to a certain degree. This 
may be proved by considering that it is well established for 
waves in the open ocean to have a maximum waveslope of about 
16 degrees. Waves much steeper than 16 degrees are unstable 
and break up by forming whitecaps or similar turbulence. 

However, rolls of 20 or 30 degress or more are common even in 
moderate seas. Thus, for the normal range of ship parameters 
and dimensions, a resonant condition occurs within the fre- 
quency range of wave exciting forces. 

To better understand the implications of being in a 



resonant condition, consider the ship as a simple second 
order system in roll, with the coordinate <j> representing 
roll angle. 

(Kp-Ix) "0 + Kp<J> + K4> = - (Wave excitations) 

The natural resonant frequency of the system is w_=^ / Kcp 

8 V ( k p- ix ) 

and it may be recognized that the amplitude of the roll motion 
at the resonant frequency is proportional to Vk^ (Kp-Ix) /Kp 
Hence the magnitude of the rolling motion at resonance will be 
most strongly dependent upon the damping parameter Kp. Some 
analyses of ship motions have assumed a linear damping para- 
meter for the ship, i.e. Kp = a constant. This is a good 
approximation for small motions or for specific conditions, 
but in the general case we expect Kp to be non-linear. Since 
roll damping involves viscous friction, we expect that there 
will be a Reynoldls number dependency. Also, because rolling 
creates waves, there is added damping due to the radiation of 
wave disturbances from the ship. This radiation of wave energy 
has been shown to be frequency dependent in a complicated 
fashion. 

Another problem is that rolling motion is not indepen- 
dent of other ship motions and therefore cannot be reasonably 
expressed as a simple second order system. For example, a 
ship at rest has been caused to assume an angle of roll which 
is different from its usual upright position. This roll angle 
causes the ship to have an unsymmetrical underwater shape, 
which induces a righting moment (assuming a stable ship) . 
Obviously, any pitching motion would change the underwater shape 
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still further, unless the ship were prefectly spherical, 
and would induce a change in the value of the righting moment. 
It can be shown by similar arguments that roll is coupled to 
all of the other ship motions and displacements. However, 
extensive studies have shown that the degree of coupling is 
negligible in all instances, except for sway rate and sway 
acceleration, assuming a "normal" surface ship. 

This paper will take an engineering approach towards 
analyzing rolling motion. The best available information on 
how to model non-linearities will be used in a workable 
system which allows calculations to be done in a reasonable 
amount of time. This approach was taken because of the 
iterative nature of the design process for Anti-Rolling tanks, 
which in turn is one step in a larger iterative design process 
for the entire ship. Therefore the analytical model of the 
Anti-Rolling tank had to be flexible enough to be "turned 
inside out" into a systhesis model by which a large range of 
parameters could be investigated by some search technique. 

This factor was the primary motivation behind taking an engin- 
eering approach to a sufficient design rather than a purely 
analytical approach towards a perfect dynamic model. It is 
hoped that this paper possesses a good enough balance between 
numerical estimation and reality to shed a little more light 
on a very complicated subject. 

2 . System of Units 

The English system of units is used throughout this 
paper and in the computer program. The following English units 
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are assumed pertinent to all applicable quantities unless 
other units are indicated. 

Time: seconds 

Distance: feet 

Mass: slugs 

Angle: radians 

2 

Gravity constant: 32.2 ft. per sec. 

Displacement: long tons (2240 lbs) 

Force: pounds 

Moment: foot-pounds 
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II. Derivation of the Equations of Motion 
1. General Discussion 

This chapter presents a detailed derivation of the 
equations of motion of the ship-tank system. The form of 
the equations is highly dependent upon the coordinate system 
chosen, which makes it important to understand section 2 and 
figures 1 and 2 before continuing with the text. For those 
who are interested only in the results, section 9 presents 
the final form of the equations and section 10 gives their 
non-dimensionalized equivalent. It is suggested that readers 
who are not interested in the derivations should still read 
section 2 before skipping to the results. This will ensure 
the reader has a reasonable feeling for the basic nomenclature 
of the equations before facing them. 

The derivation of the equations of motion begins with 
Newton's second law of motion, action equals reaction. All 
of the terms which make up the reactive forces and moments are 
first listed. These terms are each broken down into their 
dependent variables. An expression for the reactions is then 
formed by taking the first partial derivative of the partic- 
ular force or moment with respect to each of it's dependent 
variables. If the indicated partial derivatives were all 
constants the reaction terms would be exactly linear. Each of 
the partial derivatives is examined for this property, simpli- 
fications are made where possible and in several cases secondary 
parameters are introduced to further define particular coeffi- 
cients. Enough documentation of relevant assumptions is made 
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to be able to change particular coefficient values whenever 
the state of the art permits. The excitation forces and 
moments are handled in the same fashion except that they 
have independent as well as dependent variables. 

The last section of this chapter gives the summarized 
results as a set of simultaneous equations which are in non- 
dimensional form. This section is the key reference in 
understanding the rest of this paper. This is because the 
non-dimensional coefficient parameters found there constitute 
a kind of universal set which not only appears in most other 
literature on Anti-Rolling tanks (often with the same symbology) 
but also is used continuously throughout the rest of the paper. 

2 . Coordinate System 

The coordinate axes for the ship are through the center 

of gravity of the ship and lie along the principal axes of 

inertia. They are x,y, and z. The reference coordinates 

with respect to the earth are x , y , and z ; positive z 

o J o o o 

direction is in the direction of gravity. 




The tank coordinate, n , is relative to the ship and is 
positive in the direction shown on Figure 2. All of the symbols 
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shown on Figure 2 are in their positive sense. 

The intent of this figure is to show the relative 
positions of points G, B, and 0. G denotes the center of 
gravity of the ship with a filled Anti-Rolling tank (cross 
duct blocked) and serves as the reference for the ship 
oriented x-y-z axes. B, the center of buoyancy according to 
Archimede's principle, is assumed to be the point of appli- 
cation of the wave excitation forces and moments. 0 is the 
origin of the tank coordinate n . Also shown is the waveslope 
angle 4 > and some of the basic Anti-Rolling tank geometry. 
Figure A-l in appendix A gives a more complete picture of the 
Anti-Rolling tank by itself. 



ship hull 
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Figure 2. X-Z Plane, Looking Forward, Indicating Positive Direction' of Axes, Forces 



3. Reaction Forces and Moments. 



For the ship-oriented coordinates a force and moment 
balance is done. K denotes a roll moment about the x-axis, 
Y denotes the force in the y direction, and E denotes the 
tank moment about an axis parallel to the x-axis. 

Z (Moments in <J> direction) = K = I<j> 

Z (Forces in the y direction) = Y = mv 

For the tank coordinate, reference (4) uses a Lagrangi, 
description of the potential, kinetic, and dissipation 
energies associated with the tank to arrive at its equation 
of motion. For the sake of brevity, only the results, 
which are well known, will be used in this paper. 

Assuming a full tank with blocked ducts, there will 
be no interaction between the tank and roll or sway. In 
all other cases, there will be coupling between the three 
coordinates. This leads to an expression for K and Y as 
follows : 

K = K(tank) + K (hydrodynamic) + K (excitations ) 

Y = Y(tank) + Y (hydrodynamic) + Y (excitations) 

The hydrodynamic forces and moments are those imparted 
to the body of the ship by virtue of the ship motions 
occuring within and at the surface of a dense fluid (sea- 
water) in the presence of gravity. The tank forces and 
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moments are caused by the tank fluid shifting its weight 
transversely in the ship. Free surface waves at the surface 
of the ocean are the cause of excitation forces and moments. 

This separation of terms in the force and moment 
equations may be carried further by noting which variables 
correspond to tank, hydrodynamic or excitation portions. 

K or Y (tank) = K or Y (n, n, n) 

• •• 

K or Y (hydrodynamic) = KorY (<p, <f> , <p , v, v) 

K or Y (excitations) = K or Y (<f0 

The usual process of obtaining a first order model 
of a system is to take the first partial derivatives of each 
quantity (K or Y in this case) with respect to each of 
its arguments. After following the above procedure, 
combining terms where possible, and separating all the reactive 
forces and moments, including those of the tank, from the 
excitations, the roll and sway equations are obtained. 

(roll) (KvD + Kv) v + [ (K£> - lx) D 2 + KpD + K4> ] <f> + K(tank) = 

-K (excitations) 

(sway) [ (Yv - m)D + Yv]v + YpD 2 + YpD + Y<f>)<{> + Y(tank) = 

-Y (excitations) 

where: D = d/dt 

Kv = 3/ ( 3 v) K Yv = 3/ ( 3 v) Y 

Kv - 3/ ( 3 v) K etc. 

Kp = 3/(3p)K 
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Kp = 3/(3p)K 
K <|> = 3/ ( 9 <(> ) K 

The first order equation of .motion for the Anti-Rolling 
tank, as found in reference (4) is: 

(tank) pA Q R 2 l2Dv - (S"D 2 + 2g) <f> - (S'D 2 + B.pD/2 p A Q R 2 + 

2g) n] = “Eq ( t) 

S' and S" are geometry factors which are calculated 
in Appendix A. 

The coupling terms between roll, sway, and the tank are: 

K(tank) = Knn + Knn + Knn 

Y (tank) = Ynn + Ynn + Ynn 

Since the coupling terms must be symmetric in any energy 
conservative system, the tank equation gives us directly 
the indicated partial derivatives in the roll and sway 
equations . 

Kn = -p A 0 R 2 S 1 ' Yq = +2pA Q R 2 

Kn = 0 Yn = 0 

Kn = -2pgA Q R 2 Yn = 0 

Using these results and taking the case of a purely 
passive Anti-Rolling tank (E Q (t) = 0) the linear equations 
of motion for the ship- tank system are formed. The result is 



a fifth order differential equation which is linear if the 
assumption is made that the coefficients of all the terms 
are constant. 
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(roll) (KvD + Kv) v + [Kp - Ix)D 2 + KpD + K<J>](f. - pA Q R 2 (S"D 2 + 

2g)n = -K 0 (t) 

(sway) [ (Yv - m)D 2 + Yv] v + (YpD 2 + YpD) <J> + 2pA Q R 2 D 2 n = 

-Y 0 ( t) 

(tank) 2pA Q R 2 Dv - pA Q R 2 (S"D 2 + 2g) 4> - pA 0 R 2 (S'D 2 + B^/ 
pA 0 R 2 D + 2g)n = -E Q (t) = 0 

In order to add a control system to the tank (mathematically) 
it is necessary to set the right hand side of the tank 
equation, -Eq (t) , equal to -A Q Rp(t); p(t) is the pump 
pressure developed across a control effector in the tank 
duct, with a (+)p direction in the ( — ) n direction. A 
separate equation for the dynamics of the pump must be 
included in the equations of motion. A simplified form of 
p(t) which should be sufficient in the general case is to 
consider that any control effector can be characterized in 
a gross sense by an equivalent first-order time lag. This 
leads to the equation: 

T p p (t) + p (t) = p c (t) 

Pq (t) = f (<J> , $ , v, n / n f >/0 

Pcft) = fgi'i’ + 92^ + 93n + g 4 n + g 5 v + g 6 <|<, | p c | < p max ' 

Wmax' 1^1+ +•••+96^1 > ?max 
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The expression for p c (t) was chosen based on linear 
regulator theory which allows feedback of all the state 
variables, plus the excitation $ . In practice, \i> is 
random and unmeasureable (economically) , and of the state 
variables, the roll and tank variables are the easiest to 
measure and the most strongly coupled. This reduces the 
controller equation down to: 

P c (t) = g^ + g 2 ^ + 93n + g 4 n, |p c I < p max 
P c (t) = ± Pmax ' I ^l* + + ^31 + g 4 M 1 Pmax 

A useful equation which must be considered along with 
the above set of state equations is an estimation of the 
power requirements of the control effector. Using the 
relationship that power is the product of pressure differ- 
ential and volume flow, and taking into account pump efficiency, 
power requirements are calculated by: 

Pump Power = (1/ep) p (t ) A Q Rfi 

where: e - 0.5 - 0.8 

ir 

4 . Simplifications of the Hydrodynamic Coefficients . 

Assume that the roll and sway motions of the ship 
as seen at the center of buoyance are uncoupled; that is, 
a roll motion about the center of buoyancy induces no 
swaying and likewise a sway motion at the center of buoyancy 
induces no rolling there. It may be easily imagined that 
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this is nearly true. Using this basic assumption, many 
simplifications in the terms of the hydrodynamic coefficients 
are possible. For example, a roll motion about the x-axis, 
which passes through the center of gravity, will be seen at 
the center of buoyancy as a roll motion plus a sway motion. 
The hydrodynamic resistance to sway at the center of buoyancy 
is felt as a sway force, but does not couple back into roll, 
with this assumption. Using the (') notation to indicate 
that a particular coefficient is taken at the center of 
buoyancy instead of being referenced to the center of gravity 
(which is the center of the roll and sway axes) the values 
of the following hydrodynamic coefficients are determined. 



Y* v 


= Yv 


Yp = 


Y'p = 


-z h Y V 


Y ' v 


= Yv 


Yp = 


Y'p = 


-Z h Yv 


K' 4 


= K<f> 


Kv = 


K'v = 


-z h YY 






Kv = 


• 

K'v = 


-z h Yv 



Introduction of the above relationships shows which 
hydrodynamic coefficients must be known as a minimum in 
order to solve the equations of motion of the ship-tank 
system. Notice that these coefficients are all defined 
as being determined with respect to the x-y-z coordinate 
system about the center of gravity of the ship. 

K4> = Static righting moment characteristic of the ship. 

Kp = Linear damping in roll. 
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K]5 = Added inertia of the ship in roll to accomodate water 
flow. 

Yv = Linear damping in sway. 

Yv = Added inertia of the ship in sway to accomodate 
water flow. 

5. Tank Coefficients . 

The parameters which determine the coefficients of the 
tank variables in the linear equations of motion are: 

pA Q R 2 = Characteristic of tank size. 

S" = Function of tank geometry = ^(d/Rjdl 

S' = Function of tank geometry = / L (AQ/a)dl 

B.J. = Tank linear damping. 

S' and S" are the weighted moment arm and sectional area, 

respectively, of the tank fluid. They are familiar terms 

. . 2 
m literature on hydraulic engineering. pA Q R has the units 

of mass times a moment arm and is a characteristic tank 

size or inertia term. 

The above tank parameters are all free to be chosen . 
However, it can be noticed that S" and S' fall within the 
following range of possible values. 

1 < S ' < «o 

-1 £ S" < 1 

| B | £ 0.05(/2gS r pA Q R 2 ) 
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For the case of S' getting very large, the tank dynamics slow 
down due to the increased fluid inertia. The extreme case 
is the equivalent of the tank being just an additional 
weight aboard ship. The sign change of S" is accomplished 
by changing the center of gravity of the tank relative 
vertically to the center of gravity of the ship. A negative 
sign indicates a high tank, a zero value occurs near the 
center of gravity of the ship, and a positive sign is caused 
by a low tank. 

The damping parameter of the tank, B^., is free to be 
chosen as long as it exceeds the approximate minimum shown 
above, say about 5 percent of critical damping. Later 
analyses will show that a good passive tank will incorporate 
far more damping than this. Otherwise, undesirable secondary 
resonance peaks appear in the frequency response of the 
ship-tank system. Therefore, for practical purposes may 
be considered a free parameter with the stipulation that 
in the final design of the tank suitable damping devices such 
as sharp corners, baffles, etc. are used to achieve the 
desired degree of energy dissipation. B^_ may be calculated 
or predicted accurately enough with standard hydraulic 
engineering 'pipe friction' tables. 

6 . Effector Pump Parameters. 



The maximum efficient effector size can be easily 
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estimated for any particular tank geometry, since the effector 

pump is installed in the cross ducting. This will enable 

the choice of a particular sized pump which we have chosen 

to characterize by its maximum pressure, p , and lag time, 

Tp. There is no particular requirement for P max r except 

that p max = (Power, max. ) ep/A 0 Rri max. For T however there 

P 

is the requirement that the pump time lag must be less than 
the natural resonant frequency of the tank for the pump 
to have much useful effect. A pump with more lag will 
introduce instabilities into the control system. 

7 . Non-Linearities . 

The hydrodynamic coefficients obviously all depend upon 
ship speed because they are caused by changing the flow 
field of the water about the ship. Of these, the damping 
.parameters, Yv and Kp, are the most speed dependent because 
they involve the dynamic lift in the horizontal plane generated 
by the ship hull serving as an inefficient lifting body. 

An angle of attack is imparted to the mean water flow by 
the relative motion between the ship and the water in the 
horizontal plane determined by the vector addition of sway 
rate and forward speed. The resultant lift force is felt 
as a resisting force in the sway direction. The static 
righting moment coefficient, K<j> , is the least speed dependent 
because the only change from zero speed conditions is that the 
constant water pressure surfaces, such as the air-water inter- 
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face, are no longer planar, but are slightly folded into a wave 
pattern. Based on these arguments, and with no further proof, 
it was decided to adopt Kp and Yv as ship speed dependent, and 
the other hydrodynamic coefficients as speed independent. 

Because the magnitude of the roll angle at resonant fre- 
quenty is much larger than its average value, the value of 
does not remain reasonably constant, say within 5-10%, within 
its expected range. Since the well known curve of static K 
versus 0 is anti-symmetric, it was decided that it would be suf- 
ficient to adopt an odd polynomial expression for , as was 
done in reference (7). The K^jzS term is replaced by: 

Ktftf + K030 3 

The K03 coefficient may be simply evaluated by noting that 
the value of K is zero when 0 is at the maximum static equilibrium 
angle, jz5r . This fixes the value of as -K0/sz5 R 3 . This can 

be even further simplified by approximately by about 80° or 
1.4 radians. This enables the use of K^3 = -Ktf/2 if better in- 
formation is lacking. 

Part of the. energy dissipation caused by roll motion is 
due to viscous friction resisting the relative motion of the 
skin of the ship to the water. Since viscous friction is Rey- 
nolds Number dependent, it must be roll velocity dependent. This 
component of roll damping can be calculated by standard hydraulic 
engineering 'friction factor' tables to include skin friction 
and the eddy-making resistance of the bilge keels. However, 
this usage precludes us from making the odd polynomial 
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approximation because the standard tables are all of the form 
l/2p (area) (velocity ) (friction factor) . This forces the 
adoption of replacing Kpp by: 

K pP + KpjpPIp 

The absolute value sign is necessary to make the symmetric p2 
function an anti-symmetric function. Since the KjjpjPiP 1 term is 
primarily due to transverse flow meeting the tremendous flow 
obstruction of the bilge keels, and not for most cases much 
influenced by forward speed, it was decided to let this quantity 
remain independent of ship speed. 

Another form of energy dissipation by rolling or swaying 
motion is found in the formation of surface waves which radiate 
outward from the ship. This energy dissipation has been esti- 
mated, using potential flow theory, for simple ellipsoidal 
bodies which resemble the underwater form of ships, and has 
been found to depend on ship speed and the frequency of motion 
(assuming harmonic motion). As would be expected, sway gener- 
ates larger waves for an average hull form and therefore this 
effect is more pronounced for sway. Thus, it was decided to 
approximate sway damping as: 

(Y vl (to,U) + Y v2 (U) ) v; 

where Y v ^ is the component associated with wave formation and 
Y V 2 is the component associated with the hull acting as a lift- 
ing surface in the horizontal plane. Roll damping may be well 
enough approximated by: 
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KpP + Kpjpptpl; 

where Kp is calculable by strip theory, or may be estimated 
from Appendix B as about .05 times critical damping. Kpp, can 
be estimated from friction factor tables at an average Reynolds 
number. If the above seems inconsistent, it should be recalled 
that roll damping is much less than the amount of sway damping 
that is coupled into roll by the equations of motion. Thus, 
the extra care is taken in the estimation of sway damping vice 
roll damping. 

To determine the above coefficients exactly requires that 
extensive model testing be conducted. Since it is envisioned 
that the method proposed by this paper is done at a stage of 
design prior to any such detailed knowledge of the ship, the 
coefficients must be arrived at either by empirical data 
derived from previous similar designs, or else by estimation 
from their theoretical values for similar geometric forms. An 
appendix to this paper gives detailed examples of how to make 
these estimations. 

. As far as the tank non-linearities are concerned, the 
restriction on the maximum tank angle, n, by the boundaries of 
the tank is the most important of the non-linearities. The 
coupling parameters S&nd S" remain almost constant for any tank 
geometry as long as centerline symmetry is preserved. However, 
a saturation effect is observed whenever the fluid in the tank 
reaches the tank top because at that point no more fluid transfer 
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athwartships is allowed so that the tank moment reaches a maxi- 
mum. To account for saturation in the equations of motion there 
are two methods. An equivalent linear 'spring constant' may be 
used; to do this change S' such that the same energy storage 
for a given tank angle is achieved for a fictitious unsaturated 
tank as is realized in the real tank. Another method is to 
change the equations of motion at the point that the tank angle 
saturates. After saturation, consider that n = n = 0. The 
normal tank dynamics resume again when the fluid can run 'down- 
hill' • This condition is met in the static case for n + 0 > 0 . 

In the dynamic case, the dynamic head must be included along 
with the pump pressure, if present. The condition for resumption • 
of normal tank dynamics then is: 

2 p A 0 r2dv - pA 0 R 2 (S"D 2 + 2g) jzS - 2pA Q R 2 gn 
+ A 0 Rp ( t) > 0. 

The tank damping parameter, B-(-, of course has non-lineari- 
ties associated with it due to the nature of the damping devices 
used in the tank. However, as mentioned earlier, because we 
are interested in tank design, B-(- can be considered a linear 
coefficient independent of the other variables with the stipu- 
lation that the amount of linear damping must be later designed 
into the tank by some means. It is reasonable to expect that an 
amount of damping in the same order of magnitude as critical 
damping is easily achievable. 
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8 . Excitations 

In order to keep the exciting function one-dimensional, 
it was arbitrarily decided to reference all the excitations to 
the waveslope. This was mostly done to be able to compare 
results with previous work. The waveslope used was not the 
surface waveslope, but an average value which the ship senses 
(remembering that surface waves are a local phenomenon which 
decay exponentially with depth) . This average value is often 
taken at the half-draft of the ship, but to be consistent with 
hydrostatics as well as hydrodynamics, I chose to use an aver- 
age value as the value at the center of hydrostatic buoyancy, 
which is nearly the half-draft. Also, this places the ex- 
citations at the center of buoyancy, and in accordance with a 
previous assumption about the non-coupling of sway and roll at 
this location, this neatly separates roll and sway excitations. 

It is now necessary to review linear wave theory for 
simple harmonic ocean waves in deep water. A later assumption 
to be made is that the real ocean is built from a spectrum of 
these simple harmonic waves. 




particle motion 



e a Surface ele- 
I „ vation 



fu) e Horizontal wave 
a velocity 




Figure 3. Wave Profile 
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Figure 3 is arrived at from the equation of surface ele- 
vation of an harmonic water wave , given in reference 3, which 
corresponds to a surface of constant (atmospheric) pressure. 
Horizontal wave velocity, acceleration, and waveslope are found 
by taking the appropriate partial derivatives, 
e = e a cos(ks-wt) 

v w = de/dt + 90° ! =~u>e a cos(ks-oot) 
v w = 9v v/ /3t =~io 2 e a sin(ks-tot) 

= 3e/9s = -ke a sin(ks-iot) 

An important factor to notice from the above math is that 

at any given point in the wave tf» and v w are in phase with each 

other but are out of phase with the wave horizontal velocity v w 

For simplicity, the wave variables may be expressed in 

complex form, with the implication that the quantity indicated 

is the real part of the complex expression. 

e = e a exp (i (ks-wt) ) 

* . o 

v w = -we; v w = l w e ; if) = ike 
v w /ifi = -w/ik = ig/w 
v w /if) = io 2 /k = g 

The last two relationships fix v w and v w as functions of if) 
Assuming that the excitations occur at the center of 
buoyancy, and using the same prescription as before: 

K 0 (t) = + K' Vw v w + K’v w v w 

Y 0 (t) = Y^if) + Y' Vw v w + Y'v w v w 
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There is an unrealistic assumption in the above derivation, 
however, in that the horizontal wave motions are not impinging 
on a static ship. They are actually influencing a moving ship, 

which suggests that for the equation to be valid, v w and v w 

• *1 

must be replaced by v w - v' and v w - v, respectively. 

By making the assumption that the hydrodynamic coefficients 
of the water moving past the ship are the same as for the ship 
moving past the water, and by adopting the classical assumption 
of Frouctethat = K^, all of the coefficients are determined. 
Notice that Y^ vanishes due to a previous assumption, but that 
K'v w and ^ ' v w cannot vanish unless the moment arm, z^ , between 
the center of buoyancy and the center of gravity (the reference 
for the x-y-z axes) vanishes. With all of the above assumptions 
and conditions included, the excitations may now be found. 

v' = v - z hP 
V* = V - Z^p 

K 0 (t) = + zhYv(vw-v + zhp) + zhYv(vw-v + zhp) 

Y 0 (t) = Yv(vw-v + zhp) - Yv (vw - v + zhp) 

All of the terms in the excitations, except for those in- 
volving tjj,v w (i{>) and v w (<jj), may be moved to the left-hand side 
of the equations of motion when introducing this expression for 
the excitations into the equations of motion. When combined 
with like terms, this doubles the value of some of the hydro- 
dynamic coefficients, and introduces coupling terms into some 
others. Also, the excitation expression added two new parameters 
instead of the one desired — the waveslope, and the phasing 
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difference between waveslope and horizontal wave motions. 

9 . Equations of Motion 

The results of the previous sections are summarized to 
completely form the equations of motion for an anti-rolling 
tank . 

(roll) [(Kp-Ix)D 2 + KpD + - zh [YvD + Yv 1 (oj,U) + Yv 2 (U)]v 

- P A 0 R 2 (S"D 2 + 2g) n + Kp^pfp - = -K 0 (t) 

(sway) zh [YvD 2 + (Yv]_U,U) + Yv 2 (U))D]^ + [(Yv-m)D 

+ Yv 1 (w,U) + Yv2(U)]v + 2pA Q R 2 D 2 n = -Y 0 (t) 

(tank) -pA 0 R 2 (S"D 2 + 2g)^ + 2pA Q R 2 Dv 

— pA 0 R 2 (S 1 D 2 + (Bt/ pAqR 2 ) D + 2g) n = 0; or- AoRp(t) 

(effector) (TpD + l)p(t) = (g^ + g2D)0 + (g 3 + g 4 D)v, 

l z gi x ill p max ± p max< otherwise 

(excitations) K 0 (t) = -Ktfty + zh (Yv^ + YV 2 ) (v w - v + zhp) 

+ zhYv(v w - v + zhp) 

Y 0 (t) = - (Yvj_ + Yv 2 ) (v w - v + zhp) - Yv(v w - v + 

zhp) 

• •• 

where: n = n = 0and|n| = nmax, whenever | n | exceeds nmax. 

• • • 

the condition for returning to n and n free is 
whenever : 

2pAoR Dv _ pAoR 2 [s" D 2 + 2g - 2pgAoR 2 n + AoRp(t) > 0; 
with the initial condition n = + n m ax anc ^ n = 0 . 
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These equations may be linearized by dropping the terms 
involving 0 and p!p| , and the restriction on tank angle and 
pump pressure. A more convenient way to express these equations 
is to non-dimensionalize them; the results in non-dimensional 
form represent a family of solutions which help to give an in- 
sight of the general design problem. 

10 . Non-Dimensional Equations of Motion 

Because further analysis in this paper uses mostly the 
linearized version of the equations, only the linear terms are 
given below. A paragraph at the end of this section gives the 
non-linear terms as well. 

Dimensionless variables were chosen as follows: 
t = w s t 
v* = v/o> s B 
p* = P (t) 2pgR 

fi , r\ , ip in radians 



D = d ( )/dt 
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The linearized equations, with complex notation for the 
excitations, are: 



(roll) 



(sway) 



(tank) 



[(1 + *b 2 Yi)D 2 + (K s + A b 2 K y )D + 1 ]0 
- 2A b [ Y iD + K y ]v* + A t (D 2 /y st 2 + 1) n 
= [(1 - A b Ay) " iA b A 0 ] 'J'(t) 

-2 Ab [ YqD 2 + K y D]0 + (y 2 D + 2K y ) V* 

- AtD 2 / yy b 2 n = (Ay + iAo) ifi ( t ) 

At(D 2 /y st 2 + 1)0 - AtD/yy t 2 v* 

+ At(D 2 /y t 2 + KtD/yt + l)n = 0; or At p*(x) 



(effector) 



2pgR(u) S Tp + l)p* = (gi + u) S g2D)0 + (g3 + o>sg 4 D) n 



Coupling constants 
At =-2pgAoR 2 /K (Z ; 

Ab = zh/B 
Ay = YvgB/K^ 

Ao = (Yvi + Yv 2 )gB/coo K 0 

Inertia constants 
Ax = B 2 Yv/(Kp - I x ) 

A 2 = B 2 (2Yv - m )/(Kp - Ix) 

Spring constant (ship) 



Damping constants 
K s = K p /l/K0(Kp - I x ) 
K y = B 2 YvVK0(Kp - I x ) 
K t = B t /pA 0 R 2 V 2 g s ’ 

Frequency constants 
y b = ut/ws 
y s t = u st/“s 
yyt = Vg/B/ «s 



= -mgGM 
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In order to preserve the non-linear equations, the follow- 
ing terras need be added to the indicated equation, and also the 
restriction on tank angle and pump pressure applies. 

(roll) Kqppl + K3/ 3 

where: K q = (u s 2 %p,) 

K 3 = 1/2 , or 1//r 2 

There is one great idiosyncracy of these equations which 
would pass unnoticed were it not for the presence of the io 0 term 
in the wave velocity - sway coupling constant, \ Q . Here, w 0 
refers to the absolute frequency of the waves relative to the 
earth. The rest of the dynamics will occur at a steady-rate 
frequency which is the encounter frequency that the ship senses 
in the waves. This encounter frequency is a function of both 
ship speed and relative heading angle, and is further examined 
in the appendix on excitation spectra. The peculiarity is that 
the ship is sensitive to both relative and absolute excitation 
frequencies, which is a direct result of the choice of sway 



excitations . 



29 



III. Excitation Spectrum 
1 . General Discussion 

The preceeding chapter completely defined the anti- 
rolling tank problem with the exception of the single inde- 
pendent variable ip (time) . This chapter is dedicated towards 
obtaining a realistic expression for ip. Ocean waves are 
essentially a random process. To account for the statistics 
of the wave process it was necessary to adopt an accepted 
empirical spectral (implying frequency) distribution function 
which had the disadvantage of being an elevation distribution 
instead of a slope distribution. The main part of the text 
of the chapter deals with making corrections and additions 
to this distribution function. 

It is first changed into a waveslope distribution by use 
of linear gravity wave theory. Correction factors are derived 
for the amount of submergence of the ship's hull, for the 
relative heading of the ship with respect to wave progression, 
for the wave encounter frequency shift — akin to the familiar 
doppler shift — and for the shape of the ship's hull. Because 
of the highly speculative nature of some of the assumptions, 
particular attention is paid to their expected range of 
validity. All of the factors are then combined to form a 
complicated expression for the excitations which can be 
inserted into the equations of motion to complete the model 
of the ship-tank system in ocean waves. 
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This strategy was adopted because it lumped as many of 
the parameters which are independent of the anti-rolling tank 
design as possible into a single expression for excitations. 
This could then be used as an input forcing function. The 
resultant form of ip' (time) had the peculiar but advantageous 
properties of being "steady state", i.e. a combination of 
regular harmonic waves, and linear while only requiring that 
the various auxiliary component parameters be analytic. Thus 
the final derived expression for the excitation spectrum ip 
included most of the design independent parameters as well as 
the statistical distribution of waveslope. 

2 . Waveslope Distribution 

In order to describe the wave excitations of the sea 

surface it was decided to use the empirically derived Pierson 

Moskowitz wave amplitutde spectrum. This is a standard 

analytical formulation which seems to adequately describe the 

real life wave amplitude spectrum for fully developed seas. 

The statistical distribution of wave amplitude over the band 

of wave frequency (0 to °°) is assumed to depend on just one 

. . . 1/3 

parameter, significant waveheight, H . Significant wave- 

height is commonly used as a measure of sea state. To be 

2 

more specific, e /6w, wave amplitude squared per unit fre- 
quency bandwidth, is a function of frequency only for any 
given sea state. For a more complete description of this 
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spectral density function, its implications and proven uses, 
consult references 3 and 11. 

e 2 /6to = 0.0081 g 2 to 3 exp(-8/to 4 ) 

8 = 0.0324 g2/(H 1/3 ) 2 

2 

g = 32.2 ft/sec. 

Using linear wave theory, it was previously shown that: 

| ip | = | 3e/8s | = ke 
where: k = w 2 / g 

Making the appropriate substitutions, the spectral den- 
sity function for waveslope is found. 

.2/* v 2 o , v -2-1 (-3 /oj 4 ) 

\p /6 to = k (to) = ag w e 

The mean square value of the waveslope spectrum is found 
by integrating the spectral density function over the entire 
frequency range. 

2 00 2 

^ = / k S _ (w) dio 

0 
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The root mean square value (RMS) of wave slope is found 
by 4 ; rms = J ' } i )2 • By the statistics of a random process with an 
assumed Rayleigh distribution, the average of the 1/N highest 



values of ^ are given in reference 1: 



<r /LN = f (i/N) 


• l/j 

r rms 


l/N 


f (1/N) 


1/1 


1.25 


1/3 


2.00 


1/10 


2.55 


1/100 


3.34 


1/10 6 


5.27 



It should be pointed out that the above numbers refer to 
the waveslope angle measured from the horizontal, and not to 
the "crest to trough" condition usually associated with 
waveheights . 

Following the convention that "significant" means the 
average of the 1/3 highest values of the process, the signi- 
ficant waveheight is: 



= 2. Oip 



rms 



This represents an acceptable peak value, for engineer- 
ing purposes, of the random waveslope process. Notice that 
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the mean square value has a different relationship in this 
random process than for the deterministic case in which peak 
value is /2 times the RMS value. 



3 . Depth Correction 

The surface elevation of the sea is a constant pressure 
surface. This pressure decays exponentially with depth until 
at a certain depth, the surface waves are not felt. By linear 
theory, this exponential decay factor is e U / where z q 

is the depth from mean surface elevation. 

Since the center of buoyancy of the ship is the point 
through which wave excitations are induced in the static case, 
and for small amplitude excitations (i.e. nearly static case), 
a reasonable approach towards excitations is to apply them 
at the center of buoyancy. The only alternative to this is 
to carry out the integration of the pressure distribution 
over the entire immersed area of the hull. By the above 
reasoning, the spectral density function which the ship feels 
may be approximated by: 



-2tbio 2 /g , . 
e ^ (w) ; 



where t, is the draft to the center of buoyancy, 



This correction factor for depth is shown in figure 4. 
Notice that for super ships, with drafts of 80+ ft and t^ of 
about 40 ft, waves are practically not felt in roll 
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excitation. For larger ships, say 20,000+ tons, all wave 
lengths less than about 50-100 ft. are unfelt, and are non- 
interactive with the ship. 




20,000' 200' 50.6' 22.5' 



Figure 4. Depth Correction 
4 . Encounter Frequency 

Since the ship forward speed and wave advance are both 
vector quantities, relative encounter frequency of the ship 
with the waves changes with orientation. Denote heading 
angle of the ship into the waves by 3 (3 is 180° in head 
seas). The applicable relationship is: 

a) = go (1 - w U cos 3/g) 
e o o 3 
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5 . Transverse Waveslope 

The equations of motion were developed for roll and sway 
coordinates which have entirely a transverse sense with 
respect to the ship. It can be shown that for waves progres- 
sing at some other direction than transversely to the ship 
the effective waveslope in the transverse direction is: 

ifj(3) = sin (3); 3 is 90° for beam seas. 

6 . Region of Validity of Excitation Spectrum 

Combining all the corrections to the original waveslope 
spectral density function, an expanded function is formed. 

\f> 2 /6w( a), t^ , 3) = sin 2 ( 3 ) e 2t b w s (w) 

The expanded form of spectral density is plotted on 
Figure 6 for beam seas (sin (3) = 1) for the normal range of 
expected sea states. A glance at the figures shows that for 
all ships with a draft of 5 feet or greater, the desired 
range of integration of is 0.2 _< oj _< 3.0 to establish 

i^rms * The on w impose the corresponding wavelength to 

fall within 20 _< X _< 20,000 feet. 

To determine the range of 3 possible, a visualization of 
the ship in long crested waves is helpful. 
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crest 




Figure 5. Ship in Long Crested Waves. 
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Figure 6. Excitation Spectral Density in Long-Crested Waves. 
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It seems reasonable to assume that the excitation at the 
center of buyoyancy is equivalent to that for the whole ship 
as long as the value of l|cos(8)| is less than a half- 
wavelength. Another consideration is that the roll excita- 
tion, the more significant of the excitations, is propor- 
tional to K^. A look at the terms in is enlightening. 

K. = pg (Immersed Volume) GM 

4 > 

K. = (a constant) GM 
<P 

BM = ( 2nd moment of waterplane area about x-axis ) 

(Immersed Volume) 

GM = BM + KB - KG 

The above implies that K^ is proportional to the 2nd 

moment of the waterplane area. It can be shown that for a 

3 

rectangle, the 2nd moment is LB /12. Therefore, it can be 

3 

seen that K^ is proportional to (beam) . 

Since the beam is not constant along the ship's length, 
only the wider portions of the ship will contribute to roll 
excitations. It is easy to. calculate that any section with 
less than 60% of the maximum beam will have a negligible con- 
tribution. Therefore, the ship can be considered to have an 
effective length Le which should be used in the formula 
L | cos (6) | to replace the waterline length L. An empirical 
formula, which is a conservative estimate based upon the 60% 
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max. beam criterion, was derived graphically using a Taylor 
Series hullform from reference 9. 




Figure 7. Effective Shiplength. 

The region of validity for heading angle is therefore given 
by : 



Le | cos ( 3 ) | A/2 

2 

but ... A = 2irg/w 

Applying this standard to graphical display is done in 
Figure 8 below. 
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Figure 8. Region of Validity for Heading Angle. 

To summarize, it is obvious that the range of parameter 
values for w and 6 vary significantly with ship dimensions, 
ship geometry, and sea state. It is important in each case 
to investigate whether the extended equations are still valid. 

7 . Heading Angle Correction 

From the results of Section 4, it would seem desirable 
to extend the range of heading angles further by some approx- 
imate means. An easy way to do this is with the assumption 
that the static righting moment is uniformly distributed 
along the effective length of the ship. Then, the wave exci- 
tation vanishes anytime that Le | cos ( B ) | /A has an integer value 
of 1 or greater. It can also be surmised that for non- 
integer values of Le|cos(B) |/A greater than 1, the part of 
the ship's effective length between successive crests 
receives a net zero excitation due to cancellation effects. 

It is appropriate to include this effect as a correction 
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factor to the excitation spectrum. An efficiency parameter, 
c, which is a function of Le, B, and X (w) is formed by put- 
ting the above arguments into mathematical form. 



let: X = Le|cos(B)|/X 



X is 


between 


e[Le, B, X (w) ] 


0.0 


0.5 


1.0 


0.5 


1.0 


2(1 - X) 


1.0 


1.5 


2/3 (X- 1) 


1.5 


2.0 


2/3(2 - X ) 


2.0 


2.5 


2/5 (X - 2) 


2.5 


3.0 


2/5(3 - X) 


3.0 


3.5 


2/7 (X - 3) 


3.5 


4.0 


2/7 (4 - X.) 


4.0 


4.5 


-0 




Figure 9. Heading Angle Correction. 
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The excitation spectral density function now is: 



.2 2 . 2 . 0 . -2thw 2 /g _ , v 

if/ /(Suj = e sm (B) e u ^ (m) 



Since an actual ship will probably not have a constant 
beam along its effective length, will not be uniformly 
distributed. The region of net zero excitation will not be 
quite equal to one wavelength, but will change with the 
relative orientation of the passing wave. Qualitatively, 
this will tend to introduce higher harmonics into the 
encounter frequency of the roll-inducing wave energy than 
that predicted by the above model. However, in general it is 
safe to assume that the cancellation effect does occur, even 
though the higher harmonics are produced. 

Another consideration is that the final product of 
interest, roll angle, will be determined from integrating the 
spectral density function, and taking the square root of the 
result. This set of operations makes the roll angle very 
insensitive to the exact form of the heading angle correction, 

Based on the above comments, it seemed consistent with a 
conservative approach to take the correction factor as the 
line connecting the peaks of the function shown in Figure 9. 



e = 1 . 0 , X < 0.5, 



e = 1/2X, otherwise. 
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8 . Spreading Function 

A spreading function was used to convert uni-directional 
long-crested seas (implied by S^) into a more realistic 
short crested sea surface. Following reference 3, the 
spreading function f(y) is a symetric function centered on 
the primary direction of the seas. 



2 

f (y ) = 2 /tt cos (y) , 

r /2 f (y) S rr (co)dy = 
irA U 



-■jt/ 2 £ y <_ tt/2 
1 . OS^ (w) 



The effect of this function is to smooth the long 
crested seas prediction for roll response somewhat. It implies 
that there is always a beam component of wave motion influenc- 
ing the ship even when the ship is headed into or away from the 
predominant direction of the waves . 



9 . Sea State Data 

The below listed values were taken from references 8, 11 and 
12, and are useful for design purposes. The peak frequency 
of the surface wave elevation spectrum S^ is the given. 
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IV. Solution for the Equations of Motion 

1 . Introduction 

Given that the coefficients are all known, there are at 
least three practicable methods of solution available. These 
are analog, numerical integration, and complex-algebraic 
manipulation of the linearized equations. Each of these 
methods has its own advantages and drawbacks compared to the 
others. The latter method was the one finally adopted because 
it was the most flexible for design purposes and yields the 
most useful information for effort put into using it. How- 
ever, further work with the other methods may be warranted 
depending upon the particular anti-rolling tank problem being 
investigated. The details of the complex algebraic method 
are given in this chapter. A brief overview of the results 
of the investigations into the other two methods is given in 
Appendix D. 

2 . Complex Algebraic Solution 

The use of complex algebra along with the spectral den- 
sity technique allows for a closed form solution of the 
linearized equations of motion. The basic assumption involved 
in this solution bears a close resemblance to what may be 
imagined as the real situation of the response of a ship in 
waves. It should be recalled that gravity waves have an 
advance velocity which is dependent upon frequency, so that 
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the phasing between ocean waves is essentially random. The 
waves of different frequency (or wavelength) travel along at 
different velocities. If it is assumed that the ocean has a 
constant spectral density with time, then the ship's response 
to each of these frequency components has reached a steady 
state value. Random phasing of these components accounts for 
the Rayleigh type distribution of peak values of response, 
and the principal of linear superposition neatly ties the 
whole system of equations and density functions together; the 
ship's response is the sum of the steady state responses to 
each of the frequency components. 

A mention should be made of the implications of the term 
"linearity". The basic mathematical operation for each com- 
ponent case of frequency excitation and response is linear 
because the coefficients have been determined. However, the 
expressions for the so-called linear coefficients which 
involve a> e , w , X, phasing, ship dimensions and geometry, 
tank saturation, heading, sea state, etc. have no such 
restrictions. They are required only to be analytical. 

The standard method for finding the steady state solu- 
tion to a set of ordinary linear differential equations is to 
make a complex transformation, using a complex algebraic 
operator for differentiation, solve the complex equations 
algebraically, and retrieve the information of response mag- 
nitude and phase angle from the complex expression. 
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Let ^ = ip c> e 1 * JT . Then the steady state solution must be 

of the form: <b = <f> e 1 ^ IT , g = r| e 11JT , v* = v* e 1 ^ IT , where 

o o o 

★ 

^o' V V o are com P-*- ex constants. The operator 

2 2 

D = d ( )/dx may be replaced by D - = iy, D = -y . The equa- 
tions of motion now become for a passive tank: 
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N (iy ) 

determinant [C ( iy ) ] 



The numerator, N (iy ) , is found by replacing the first 

column of the C matrix with the coefficient matrix for \b . 

Y o 

Notice that the C matrix terms will not be symmetric because 
of the choice of v* as a variable. If y* were chosen instead 
of v*, the symmetric condition for an energy conservative 
system would have been met. 

Since p = iy<{> by linear theory, the roll velocity 
response may be obtained by y|^-|, and likewise the roll 



48 



acceleration is P 2 !^!* Tank angle, velocity and accelera- 
tion are similarly determined. 

The final step in obtaining roll response is to multiply 
the transfer function <p 2 /\p 2 (or 2 /4 > 2 » or <j>* 2 /i> 2 ) by the 
excitation spectral density function ^ 2 /<5w and integrate the 
result over the entire frequency band. The spreading func- 
tion f(y) may be inserted in this step depending upon whether 
long-crested or short-crested sea result are desired. 



rms 
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V. Anti-Rolling Tank Design 
1 . Introduction 

The complex algebraic form of solution of the non- 

dimensional equations of motion was written in computer 

language to allow the numerical work to be done by a digital 

computer. Details of this computer program are given in 

Appendix C. Essentially there are two separate programs — 

1/3 

MAIN 1 which maps <p versus y^ and for one case of head- 
ing angle in long-crested seas and MAIN 2 which calculates 
short-crested sea results for any given set of input para- 
meters. In the interest of numerical efficiency most of the 
iterative steps of the design process can be done by MAIN 1, 
looking at the largest component of roll in long-crested seas 
which makes up the center frequency of short-crested seas 
results. MAIN 2 can be used later in the final stages to 
give short-crested sea roll angle predictions based upon the 
results of using MAIN 1. 

The design method basically entails "bouncing" back and 
forth between MAIN 1 and MAIN 2 looking for optimum para- 
meters of the non-dimensional equations. The steps are man- 
ually done, meaning that each step represents a separately 
submitted computer run. About 2 to 4 computer runs are 
required for a final design for each different anti-rolling 



tank . 
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The passive tank design method in section 2 suggests how 
to specify various parameters of the equations and presents a 
logical sequence to arrive at an approximate tank size. 

There is also a suggested iterative sequence to be used as a 
design aid. One of the fortunate results of the subsequent 
chapter on design examples is the collection of equations 
shown on Figure 12. It is recommended that this worksheet be 
kept as a tally sheet and guide when performing the iterative 
steps . 

The active tank design method given in section 3 is 
mostly a discussion of possible changes to the passive tank 
design which are compatible with the computer program. Mak- 
ing these changes will result in a prediction of the perfor- 
mance of a feasible activated anti-rolling tank. 

2 . Passive Tanks 

The response of the ship-passive tank system is modeled 
by the non-dimensional equations developed previously and 
setting E Q (t) = 0. 

Tank location may be specified or a small number of dif- 
ferent locations may be investigated. The reason for limit- 
ing tank location to the absolute minimum is entirely practi- 
cal. The vertical space of the ship is divided by structure 
into a number of decks and for access and strength considera- 
tions the cross duct of the U tube tank should run across (or 
under) a deck. Also there may only be a single or a few 
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regions within the ship large enough to accept the tanks. As 
a consequence of the non-dimensional equations, there is a 
curious coupling between tank vertical placement, size, and 
the coefficients of the equation's. Since the ship's vertical 
center of gravity is shifted with any vertical shift of the 
tank, the ship's GM is reduced whenever the tank is moved 
higher in the ship. Since is proportional to GM, and w s 
is proportional the square root of , and because all of the 
non-dimensional coefficients depend on or , each verti- 
cal tank location must be treated as a separate problem. The 
design method which follows is for each separate case of ver- 
tical tank location. 

Once tank location is specified, the parameters which 

are still free to be chosen are: A. , y . , y . , K.. We know 

t t st t 

that has a minimum value of about -0.05, and that y will 
have a value close to 1.0 because the vibration absorber 
effect of the tank is maximized by tuning the tank to 
y fc = = 1*0. These four dimensionless parameters are 

dependent upon the real ship and tank parameters as follows: 

A t = A t (GM) 'P\> r2 ^ 

y t = [w g (K^) , R,S'(A o /A 1( L o /R)l 

y st = ^st [w s (K )^S"(Z t /R,L o /R)] 

2 

K-t = K t [pA Q R , geometry of tank, 
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If the ship changes its loading condition during 
expected operations such that its displacement is appreciably 
changed, then all four parameters have a A dependency, also. 

Another important parameter for design is ip max. Tank 
capacity ip max is defined as the maximum heel angle that the 
tank can impart to the ship by maximum movement of the tank 
fluid off to one side. This also represents the maximum 
waveslope angle which the tank is capable of stabilizing. To 
form an opinion of a required ip max, use must be made of 
empirical formulas. Hagen mentions the use of: 

fmax = 0- 36/Log 10 (A) radians 

A formula suggested by center of buoyancy correction to the 
spectral density function which takes the equivalent of a 6° 
surface waveslope at the center of buoyancy for the ship's 
natural resonant frequency is : 

^ = 0.10 e " (t b w s 2 /g) radians 

max 

The above formulas give an idea of the approximate capacity 
which a tank must have in order to be effective, and serve as 
a feasibility check before proceeding further with a particu- 
lar design. 

2 

The parameter = (-2pgA Q R /-mgGM) has certain char- 
acteristics which allow it to be specified by the designer. 
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2 

Rearranging the formula to read A = [2A Q R /(displaced volume 
of ship)]/(GM), notice that the numerator of this ratio is 
the standard free surface correction for a virtual rise in 
the ship's center of gravity due to the fluid being unconfined 
in the n direction. Therefore, A^ is the percentage loss in 
static stability of the ship due to the tank fluid free sur- 
face. There will be a maximum allowable A^ which ensures 
safe operation of the ship under all loading conditions, 
including damaged conditions. Since a larger A^_ is associ- 
ated with more tank stabilizing effectiveness (t|> -A^L /R) , 

3 max to 

the choice of A fc is fixed at the maximum value which allows 
safe operation of the ship. 

An important conclusion to reach from the above observa- 
tions is that the ship hull form should be designed with the 
notion that an anti-rolling tank is to be installed. This 
means that space, weight, and static stability will be such 
that the ship will be comfortable with a tank of the capacity 
indicated by the empirical formulas. The alternative choice 
of trying to fit an anti-rolling tank into an existing ship 
design which may not have all the necessary reserve space, 
weight, or static stability most likely will result in an 
ineffective stabilizing system. 

Now that A^ is known, and \ 1 ) has been determined 
t r max 

approximately, an assumption for a y value is all that is 
required to calculate tank dimensions (A q , R, L q , z^_) the 
tank parameter y st (L Q /R, z /R) , and approximate tank weight 



54 



[2pA Q R(L 0 /r + A^/A ) ] . A first guess at y fc is made by using 
Chadwick & Klotter ’ s equation, based upon the linear equations, 
neglecting sway, from reference 4. 

M t = a - A t u/u st 2 - i/p st 4 >r 1/2 

It is now possible to calculate all the tank dimensions 
and all the coefficients of the equations of motion with suf- 
ficient accuracy to consider them constant, except for duct 
area A^ and tank parameters y and K t> 

The remainder of the design problem may be stated as 
follows: Find the best combination of y and K which mini- 
mizes the rolling motion predicted by the non-dimensional 
equations of motion. This two parameter search problem may 
be solved by noticing that y must be about 1.0, and that two 
digit accuracy is the best to be expected from the equations 
of motion. Therefore, instead of a sophisticated search 
technique it should suffice to map roll angle versus y^ and 
K^_. This may be. further simplified by sub-optimizing each 

y t value for K^. This approach is justified by gentle cur- 
l/3 

vature of the function <J) (y t ,K^_) , in the region of interest, 

as found in all sample problems tried. The remainder of this 
section will investigate numerical methods to deal with this 
search problem. 

The basic strategy for solving this problem is to first 
fix upon some criterion of roll performance. Based upon the sea 
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state data already given, the author decided to minimize roll 
angle at the worst heading at operational speed in sea state 
6 (15 feet significant wave height) was a reasonable criter- 
ion. Other criteria could include roll velocity or roll 

acceleration for other environmental conditions. 

1/3 

In order to compute 4> for given conditions, a long 

1/3 

sequence of operations is accomplished--^) is twice the 

2 2 2 

square root of the integral of <p /y times <*) /<$w. This 

sequence makes vivid the reason why the complex algebraic 

method of finding <f >/\p is superior to numerical integration. 

The latter would require another integration of many steps, a 

very expensive approach since the results must in turn be 

iterated many times for a good design. 

The above discussion assumed that long crested seas were 

involved. For short crested seas another integration, this 

time over every possible heading angle, must be performed for 

the spreading function. For this reason, long crested seas 

1/3 

are assumed for mapping <{> [y^, K (best) ] . Short crested 

seas are only looked at later as an output. 

1/3 

In mapping <J> [y^_, K^fbest)] , it was found that the 

function displays very gradual curvature as was to be 
expected for an integral function. Therefore, only a few 
points are investigated during each computer run, because 
interpolation is predictable. A second order Legendre 
polynomial curve is used in an interpolation scheme to find 
a local minimum or maximum. Manual inspection of the computer 
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results will show whether further mapping is necessary. 

Warning: the interpolation function likes to extrapolate to 

find minimums or maximums equally well, but such results are 

invalid by the statement of the problem. 

1/3 

Once < p [y^, K t-(kest)] has been mapped and values of 
and for a minimum roll angle have been decided, short- 
crested sea results may be obtained. The problem then shifts 
to chosing which of the alternative tank locations is the 
best. Economic as well as performance factors will enter 
into the decision of this question. For example, a tank 
which is in a cargo hold, presumably valuable real estate, 
would probably not be chosen over one which was in a less 
costly part of the ship, such as atop an outside deck, were 
the two tanks to perform nearly as well. 

A detailed listing of the author's steps in passive tank 
design are reported as an aid in improving upon them. A 
later section will include a completely worked example of the 
method. 



I. 



Problem setup 

a) Assume a range of a. small number of vertical tank 
locations . 

b) Calculate ^ max requirement. 

c) Specify a X 

d) Assume an initial y value from Chadwick and Klotter ' s 
equation . 
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e) Determine Lo, A q , R, z , and approximate . 

f) Design the tank except for A^ . 

g) Calculate the hydrodynamic coefficients and other 
computer inputs. 

II . First Iteration 

a) Run MAIN 2 program; use the long crested wave results 
to estimate the worst speed and heading angle com- 
bination. Use y assumed, K = -0.8 

1/3 

b) Map <J> [y^, K -(-(k est )] using MAIN 1. Use y assumed, 

K initial = -0.4, Ay t = 0.5, AK^ = -1.0. 

1/3 

c) Interpolate to find <J> minimum. Check computer 
results to ensure that a true minimum occurred. 

III . Second Iteration 

a) Run MAIN 2 again to see if the worst speed and head- 
ing angle combination has changed. Also, roll per- 
formance may have degenerated at other headings and 
speeds significantly. 

b) Continue iterations if change in speed or heading 
angle warrant it. 

IV. Alternative Tank Choice 

a) Based upon the operator's or owner's criteria, select 
the best tank location. 
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V. Verification 

As a final step, short crested sea results, MAIN 2, 
should be run for the entire expected range of speed and 
sea conditions. This is to make sure that the solution 
didn't fix too shortsightedly upon some assumed worst 
condition and lessen it without overlooking its negative 
effects on some peculiar or unforseen combination. 

Another verification which should be performed during 
this and all other steps as well is to check the range of 
parameters, the range of the solution variables to ensure 
that their coefficients remain linear (or else adjust the 
coefficients to their equivalent linear values) and that 
the and K values being predicted are realizeable. 

3 . Activitated Anti-Rolling Tanks 

An activated anti-rolling tank is one which has an input 
power source in addition to ship induced excitations. A 
pressure source is used as the type of power input in this 
discussion, as opposed to a flow source. It is desired to 
regulate the pressure source in such a way as to reduce the 
rolling motion of the ship. 

The choice of a pressure source instead of a flow source 
was not entirely arbitrary. A flow source would have to be 
some sort of a positive displacement pump capable of handling 
large volume rates of fluid transfer, approximately 500 tons 
per minute as an order of magnitude. The required pump size 
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for this volumetric capacity would make it economically 
unfeasible. Another requirement for the power source that is 
imposed by economics and amply illustrated in other papers, 
e.g. reference 18 , is that it operate at a nearly steady 
power level. The peak power output for a pressure source, as 
found by standard control theory methods using linear feed- 
back, is too high to be practical. There is too much wasted 
power capacity which is only used at peak power demand. How- 
ever, the pressure pump must be built large enough to be able 
to supply the predicted peak demand if the linear equations 
are to be valid. Once further economic factor to be con- 
sidered is the competition, the alternative methods of roll 
damping. With regards to these, the gyroscopically controlled 
anti-rolling fin installation is the most serious candidate. 

In practice, a standard, well designed anti-rolling fin 
installation has been found to consume about one or two per- 
cent of the installed hoursepower of a ship, including pro- 
pulsion, due to the added effects of the fin area increasing 
the ship's total hull resistance to passage through the water 
and the power needed by the fin positioning motors. Thus, 
economics and practical engineering considerations have dic- 
tated a rather narrow range of possible control effectors for 
anti-rolling tanks — a pressure type source of power of nearly 
constant output which has a power rating of about one percent 
of the ship's installed power. 
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A further consideration of the power source is reversi- 
bility. A level power output does not imply that the direc- 
tion of pressure or flow will be constant. The pressure 
supply will obviously have to reverse itself as the ship rolls 
back and forth. This aspect of the pressure source, quanti- 
fied by a characteristic time lag to change direction, will 
have to be reasonably estimated and included in the analysis, 
and must also be reasonable in the sense that a smaller time 
lag probably represents a more expensive piece of equipment. 

It is helpful to first look at the orders of magnitude 
of the parameters needed to describe a control effector. 
Consider a ship of 20,000 tons displacement with a speed of 
20 knots. Assuming an overall lift-to-drag ratio of 150, the 
installed horsepower would be about 18,400 Hp. This would 
indicate from the preceeding arguments that an economically 
feasible anti-rolling tank can be no greater than 184 Hp. 
Compare this to the energy flow of the ship in its roll 
coordinate, assuming a GM of 6' , a 5° roll, and a roll period 
of 12 seconds. 

20,000 tons (2240 lb/ton) (6 ') (sin 5°) (1/2)- 
(1/12 sec) (Hp-sec/550 ft-lb) = 1775 Hp 

This simple estimate is enough to show that we cannot 
expect the control effector to be able to counter the effects 
of wave excitations. However, it may be used to counter the 
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resonance phenomenon associated with rolling frequencies near 
the natural resonance peaks of the ship-tank system. 

Since the passive anti-rolling tank is capable of devel- 
oping nearly as much energy flow as the ship, until the 
saturation angle of the tank is reached, it seems justifiable 
to further simplify the control effector equation to 

T p p(t) + p(t) = ±P c 

P c denotes the steady pump pressure developed across the tank 

duct, and Tp is the reversing lag time for the pump. The 

• • 

equation P c = f (<J> , <j> , h , n) implies now that the sign for P c is 

chosen according to some criterion contained in the state 

• • 

variables <t> , 4> , n , and h • Consider the same ship with an 
anti-rolling tank of 200 tons weight of fluid. The average 
fluid volume flow would be: 

0.707(200 tons) (35 ft 3 /ton H 2 0) (1/12 sec) = 142 ft 3 /sec 

For an average 184 Hp, the mean pressure is 184 Hp (550 ft-lb/ 
Hp-sec) (1/142 ft 3 /sec) = 245 lb/ft 3 =1.7 psi. This pressure 
changes the fluid level between the tanks by 3.8 ft. Assum- 
ing a 50 percent pump efficiency, the pressure differential 
would be 0.85 psi with a fluid level difference of 1.9 ft. 

This corresponds to a tank angle of about two degrees. The 
static ship roll angle for still water and constant pump 
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pressure would be A times two degrees, amounting to only a 
fraction of a degree. Thus we can now state that the addi- 
tion of such a control effector would have a qualitatively 
negligible effect on the equations of motion, except near the 
resonant frequencies. This indicates that some form of 
active damping is the goal of the control effector switching 
criterion . 

The problem of adding a control effector now becomes one 
of how to add the active damping to the linearized equations 
of motion. The tank equation is the only one of the set 

which contains p(t) for an active tank. Substituting P c = 

• • 

g^<J> + g^<p + g^n + g^r) into the tank equation, neglecting for 

the moment T (consider T small) , and combining terms in the 
P P 

linearized, non-dimensional equations of motion, the struc- 
ture of the active tank equations is revealed. 
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We would expect that the effects of g^ and g^ to be 
relatively small even if we took them at their maximum pos- 
sible values ( to develop 184 Hp) because of the small static 
values of tank and roll angle the pump is able to develop. 
Hence, making g^ = g^ = 0 is approximately correct for any 
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economical tank. The choice of g^ anc ^ g^ to use as a switch- 
ing criterion is obvious. The tank already has linear damp- 
ing which may be varied as required by adding valves or 
baffles. Negative tank damping 'is undesirable because it 
would increase roll response at resonance frequencies. The 
computer program results from MAIN 1 confirms this hypothesis. 
Adding a damping term to the ship roll-tank coupling terms, 
however, would have a desirable effect at any frequency. An 
additional tank moment caused by the pressure differential 
would be induced opposing the roll velocity of the ship. 

Thus, the controller effects may be included in the linear- 
ized equations of motion for the type of controller described 
by substituting this reduced tank equation into the non- 
dimensional equations of motion. 
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The value of g^ is determined by power considerations. 
The energy dissipation, the product of assumed pump pressure 
and tank volume flow, can be calculated and compared to the 
desired pump power in an iterative fashion until the two 
balance. The effect of pump lag time may be included by 
deducting the fraction of the period of roll that was spent 



64 



in changing direction. This factor is about 2 tt/iu - T^ for 
Tp less than 2tt/w and zero for T^ greater than 2 tt/ui . 

The above active tank equations are only useful for pre- 
dicting ship motions in a seaway which overpowers the passive 
aspects of the anti-rolling tank. Fortunately, this is the 
case when an active tank would most likely be used. The 
additional damping term contained in and the need to cal- 
culate mean pump power are small modifications to the com- 
puter program for the equations of motion. Otherwise, the 
same design method as that for the case of passive anti- 
rolling may be used. 

The type of hardware which would give the desired pump 
characteristics is a waterjet type pump situated in the tank 
cross duct. The jet velocity would be much higher than the 
mean fluid velocity, giving rise to a steady pressure inde- 
pendent of average tank fluid velocity. The presure rise is 
caused by the cross duct acting as an inefficient diffuser 
and the pressure reverse may be accomplished by switching the 
waterjet outlet from a nozzle in one direction of flow to an 
opposite nozzle. The actual pump motor would then not have 
to reverse at all, and the lag time would be associated with 
the time required to switch nozzles. 
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VI. Design Examples 

In order to validate the computer program, verify the design 
steps postulated, and investigate the possible range of the solu- 
tions, it was decided to design a series of antirolling tanks 
for a sample family of ships. The ships were chosen to be able 
to illustrate as many aspects and implications of the equations 
of motion as possible within time and effort constraints. 

Three ships were finally settled upon to show basically the 
effects of ship size, degree of sway coupling, and tank size. 

They were a 3600-ton small naval escort type, a 14,500-ton medium 
sized naval vessel which resembles a small aircraft carrier 
or amphibious landing ship, and a 20,400-ton commercial stores 
carrier. The three ships were assumed to be existing designs 
to which an antirolling passive tank was to be added. The two 
naval vessels were possessed of a large degree of sway coupling 
introduced by the parameter Xb = BG/B because they were ships 
with high centers of gravity relative to their other dimensions. 
The commerical stores carrier had only about half as much sway 
coupling, being more typical of a cargo vessel which does not 
possess the high topside payload weights of a naval vessel. 

1. Ship Size 

The principal dimensions of the three vessels were as follows 

Ship No. 1 Ship No. 2 Ship No. 3 
A 3600 14,500 20,400 



L 



408 



585 



500 
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Figure 10A Midship Section Ship No. 1 and 3, With Tank Locations 
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Figure 10B Midship Section Ship No. 2, With Tank 
Locations 



Ship No. 1 
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B 

T 

Cb 



45 

14.5 

0.4733 



Ship No. 2 
76.67 
21.33 
0.5304 



Ship No. 3 
71.4 
28.6 
0.70 



2 . Tank Size 

Based upon the empirical formula of Section V.2, an approxi- 
mate required value for tank capacity, Tinax., was determined. 

The vertical layout of the ships, by decks, is shown in Figure 10- 
Tank dimensions and parameters were calculated using the collec- 
tion of equations shown on figure 12 by an iterative process. 

The A^ value for all tanks was maximized in the following way: 

An assumed A t was used for the low tank and applied to the value 
of metacentric height, GM, for the low tank. The residual stabi- 
lity of the reduced metacentric height given by GM (1 - A ) 
was kept as a required value to be achieved by the other tank 
locations. Since the higher tank locations resulted in a higher 
center of gravity for the ship and therefore lower GM, their 
Aj. values had to be necessarily lower in order to achieve the 
required residual stability GM (1 - A^) . 





Ship No . 1 


Ship No. 2 


Ship No. 3 


¥ max. = .36/logio A 


5.8° 


5.0° 


4.8° 


¥ max. = .le“ tbw s //g 


5.4° 


5.4° 


5.2° 


GM (1 - A t) 


2.69 


5.17 


5.25 



Of the three deck locations possible for tanks on Ship No. 1, 
it was found that the tank placed on the main deck had too small 



69 



a capacity ¥ max. and ship- tank coupling A to be practical. 

Only the lower tank and the center tank were further investigated. 

On Ship No. 2 there were only two tank locations possible, 
assuming that keeping the center hangar deck area clear was 
imperative. The higher tank was slightly cramped in order to 
fit it between decks, while the lower tank should probably be 
expanded to meet the hangar deck bottom to eliminate the small 
wasted space between the tank top and the hangar deck. 

On Ship No. 3, the presence of vertical cargo access hatches 
and cargo storage limited again the tank location to two possi- 
bilities: a low tank with the cross duct underneath the cargo 
floor, or to a second deck high tank of small dimensions to 
accommodate the cross duct between cargo hatches. 

For the naval vessels. Ships No. 1 and 2, the tank fluid 
was assumed to be fuel oil such as diesel oil or JP-5 ( P =1.69 

3 

slugs/ft. ) because such vessels carry a generous reserve fuel 
supply for extended range purposes. In effect, this made the 
passive antirolling tanks a reserve fuel oil tank which is pos- 
sible to ballast with seawater in an emergency. For the com- 
mercial vessel. Ship No. 3, the tank weight was too large to 
consider it to be carrying expensive fuel oil, so seawater ( P = 
1.99 slugs/ft. ) was assumed to be the tank fluid. This re- 
duced the tank volume somewhat, compared with a fuel filled 
tank, but the tank weights of between 1-3 percent of the 
ship's total displacement was a factor against their being a 
practical commercial installation. 
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3. Operational Criteria 

The choice of speed and sea state was also made to illus- 
trate the range of solutions of the equations of motion. The 

1/3 

smallest ship was examined in a sea state 5 (H =10 ft.) 

at its operational speed of 25 knots, a typical requirement 

for an escort ship of its class. Ship No. 2 was looked at 

1/3 

in sea state 6 (H / =15 ft.) at a lower speed of 20 knots 

because it can be expected to behave better than its lighter 

counterpart, and has a lower designed speed. Ship No. 3 was 

examined at its economical cruising speed of 16 knots in a 

1/3 

storm of sea state 7 (20 ft. H ) . This is because the tank 
should be able to mitigate motion effects on the cargo in ex- 
treme conditions. 

4 . Initial Program Values 

An initial estimate of the worst heading angle, the point 
for optimizing roll angle versus y and K t , was made by taking 
the point where the peak of the excitation spectral frequency 
distribution corresponded to the higher frequency peak of the 
roll response, assumed about 130 percent of w . Later work 
showed that this was not consistently a very good first guess, 
and that a better first estimate could be gotten by running 
Main 2 program with Kt = -0.8 and taking the worst heading angle 
directly from the results. 

The initial values for heading angle and other parameters 
for the first iteration were as follows for each of the ship- 
tank combinations. 





Ship 


No. 1 


Ship 


No. 2 
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Ship No. 3 




low 


center 


low 


high 


low 


high 


BG 


9.75 


9.93 


16.5 


17.3 


8.37 


8.69 


GM 


3.36 


3.18 


7.39 


6.67 


7.00 


6.68 


fc b 


5.67 


5.67 


8.24 


8.24 


13.43 


13.43 


At 


.20 


.155 


.30 


.22 


.25 


0.075 


yt 


1.01 


1.01 


1.02 


.97 


1.01 


1.00 


2 

yst 


4.11 


4.10 


4.72 


-24.61 


4.25 


-88.3 


Kp 


-.221 


-.221 


-2.92 


-2.92 


-3.17 


-3.17 ( 


1 Y^/m | 


.67 


.67 


.53 


.53 


.88 


.88 


1 Kp/IX I 


.26 


.26 


.41 


.41 


.15 


.15 


D1 


-8250 


-8250 


-30200 


-30200 


0 


0 


D2 


-21450 


-21450 


-78500 


-78500 


-96790 


-96790 


D3 


-10300 


-10300 


-37700 


-37700 


-26800 


-26800 


FI 


0 


0 


0 


0 


0 


0 


F2 


. 684 


.684 


.684 


. 684 


.735 


.735 


F3 


1.05 


1.05 


1.05 


1.05 


1.13 


1.13 


3 


9 0° 


9 0° 


118° 


118° 


50° 


50° 


U 


43 


43 


34 


34 


27 


27 


H i/ 3 


10 


10 


15 


15 


20 


20 


KG 


18.58 


18.76 


30.4 


31.1 


23.54 


23.86 


K<J) 


-2.71 


-2.57 


-24.0 


-21.66 


-31.99 


-30.52 (X 


w s 


. 5413 


.5271 


.4456 


.4233 


.505 


.493 


Ao 


124.5 


91.3 


607 


357 


611 


214 


Ao/A ± 


5.0 


5.26 


4.61 


4.85 


3.4 


3.83 


R 


20 


20 


33 


35 


32 


32 


Lo/R 


.498 


.538 


. 303 


.286 


.349 


.313 
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Figure 11A First Iteration Tank Optimization 



High Tank 
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Figure 11B Second Iteration Tank Optimization 
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Figure 11C Third Iteration Tank Optimization Figure 11D First Iteration Tank 

Optimization 
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Figure 12. Tank Sizing Equations - Worksheet 
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Ship 


No. 1 


Ship 


No. 2 


Ship 


No. 3 




low 


center 


low 


high 


low 


high 


Weight 


84 


64.6 


506 


299 


611 


163 


Z t/ R 


.343 


-.339 


.436 


-.78 


.23 


-.67 


“st 


1.097 


1.067 


.968 


i2 . 1 


1.041 


i4 . 75 



5. Summary of Results 

1/3 

The results of calculations for <{> as a function of y t at 
(best) are shown on figure 13 for the first three iterative 
steps of looking at an assumed point in the sea state - heading 
angle space. A discussion of results for each of the three 
ships follows this. It should be repeated that the intent was 
not to design a particular antirolling tank or to apply any 
arbitrary criteria, but instead to get as broad a look as pos- 
sible at the range and implications of the equations of motion, 
and the solution method. 

Because the center tank of Ship No. 1 did not seem to offer 

any better performance than the low tank, it was not carried 

beyond the first iteration. This low tank in Ship No. 1 became 

the most difficult of the tanks to design because the plot of 
1/3 

<J> versus y t at (best) exhibits more curvature and other 
surprises than the other tanks. At each iteration, the sub- 
optimized parameters which worked well at that particular head- 
ing angle were poor at different headings. Finally, y^ was 
chosen as 1.01 by inspection of all the data, and K t was deter- 
mined by running Main 2 program for short crested seas through 
K t = -0.4, -0.8, and -1.0. The combination y^ = 1.01, = -0.8, 
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Figure 13 Ship No. 1 Low Tank 
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Figure 14A Ship No. 
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Figure 14B Ship No. 2 Low Tank 
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Figure 14C Ship No. 2 High Tank 
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Figure 15 Ship No 3 Low Tank 
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which was the first iteration guess, was just about optimum. 
Overall performance of the antirolling tank compared to the 
ship without a tank was not spectacular, but because the tank 
uses reserve fuel oil as a working fluid and is partially sit- 
uated in normal tankage space, it was a reasonably successful 
design. 

On Ship No. 2 the low and high tanks were both carried 

through to a final iterative step. The initial guess of 118° 

as a worst heading angle delayed the solution somewhat because 

1/3 

at that heading the characteristic of <f> ' did not flatten out 
as expected and also 118° was nowhere close to the real worst 
heading. Both high and low tanks performed well and it is 
interesting to note that in spite of the much smaller size, 
weight, and of the high tank it performed as well as the low 
tank. This was probably because the parameter y 2 for the 
high tank had a negative sign, indicating that roll acceleration 
pushed the tank fluid in the direction to give a positive tank 
righting moment. For h st 2 with a positive sign, the case for 
any tank below the ship center of gravity, the opposite reac- 
tion of tank fluid to roll acceleration occurs. 

On Ship No. 3 only the low tank was investigated fully, 
but it can be noticed at the 50° heading angle that the high 
tank's performance is nearly as good as that of the larger low 
tank. Again, this was due to the negative sign on y 2 . 

A relative comparison of the predicted solutions of the most 
important tank design parameters is in order. 
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Ship No. 1 


Ship No. 2 


Ship No. 2 


Ship No. 




low tank 


low tank 


high tank 


low tank 


¥ max 


5.7° 


5.2° 


3.2° 


5.0° 


X t 


0.20 


0.30 


0.19 


0.25 


U t 


1.01 


1.02 


1.20 


1.01 


V Bt 


2.03 


2.72 


i5 . 30 


2.06 


Weight/IOOA 


2.33 


3.49 


1.57 


3.00 


K t 


-0.8 


-1.0 


-1.0 


-.4 


e 


0.73 


0.26 


0.25 


0.26 



A performance parameter e (e max (tank) / max (no tank)) 

was added at the bottom of the list, but it is difficult to attach 
much relative meaning to the parameter because all of the ships 
had different operational criteria. Even if the criteria were 
the same, the size differential in the ship's displacements would 
have made comparisons meaningless except in the case of the tanks 
on the same ship. The tank capacity, 'i'max, did not enter into 
the optimization because all of the significant tank angles cal- 
culated were less than their saturation angles, the arctan (L o /R) . 

Of the parameters given above, the first four, ’Frnax, y t , and 

y s j. are the most fundamentally important because they determine 
the antirolling tank size, dimensions, weight, placement in the 
ship vertically, and they cannot be changed once the installation 
is completed. however may be readily adjusted to suit the 

particular ship and the other two parameters relate to performance 
or economic matters. 



84 



6 . Analysis of Results 

a. Tank Capacity . The question of how well the tank capa- 
city was utilized may be answered by inspection of the results 

of significant tank angle at the worst heading in short-crested 

1/3 

seas. The predicted tank angle, g ' , may be compared to the 

tank saturation angle, arctan (Lo/R) , to give a measure of how 

much capacity was used at the design conditions. Two tank angle 

1/3 

statistics were looked at. n is a good engineering estimate 

of the point at which saturation effects must be taken into 

account in calculations. ^1/100 was arbitrarily taken as an 

estimate of when tank saturation effects might just start to be 

noticed by the crew of the vessel because of the slapping of the 

1/3 

tank fluid against the top of the tank. n is an output of 

the MAIN 2 program and i s f OU nd by the statistics of the 

1/3 

random process to be 3.34/2.00 times g . 





Ship No. 1 


Ship No . 2 


Ship No. 2 


Ship No. 3 




low tank 


low tank 


high tank 


low tank 


/ /3 


9.81° 


7.07° 


7.24° 


13.48° 


1, 1/100 


16.38° 


11.81° 


12.09° 


22.51° 


arctan (Lo/R) 


26.47° 


16.86° 


15.96° 


19.24° 


g'^'Varctan (Lo/R) 


0.47 


0.53 


0.58 


0.89 


gl/100/arctan (Lo/R) 0.62 


0.70 


0.76 


1.17 


An immediate 


observation 


is that the 


low tank on 


Ship No. 3 



is the only one which is making good use of its tank capacity at 
design conditions. The other designs are all wasting space by 
building too large a tank for the requirements stated. 
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b. X . The ship-tank coupling parameter did not by itself 
seem to influence results much except as it contributed to the 
tank capacity. Ship size and tank placement had a more definite 
effect. 

c. y t . The tank tuning parameter had a secondary effect 
on the results because a higher y t corresponding to a higher io t 
is realized by decreasing the area of the cross duct, Al. The 
lighter cross duct weight allowed for more of the fluid to be 
in the wing tanks. Consequently, a change in tank tuning to a 
higher tank frequency made for a tank of slightly greater capacity. 

d. y . The ship- tank decoupling frequency parameter was 
important for the reasons previously discussed concerning it. 

A more negative value of l/y st 2 , which is accomplished by moving 
the tank higher vertically, is always associated with better tank 
performance for a tank of a given X t and capacity Tmax. 

e. Weight Fraction . This parameter can be directly related 
to the cost of the tank. A complementary parameter, the space 
fraction, was not given because the data was unavailable for 
useful ship volume. 

f. I< t . The tank damping parameter falls within the narrow 
range of about -0.4 to -1.4 for all successful designs. This 
statement is based on analysis of much more data than is presented 
in the results section. A convincing reason for the truth of this 
statement can be found by considering the critical damping for 

a single degree of freedom system. For a single degree of freedom 
system critical damping implies that there is no overshoot and 
no oscillatory response to a step input. Such critical damping 
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has the K t value of -2.0. Therefore the ratio of tank damping 
to critical damping falling in the range of 0.2 to 0.7 is rea- 
sonable if the tank is to act as a dynamic vibration absorber. A 
K t prediction of greater than -2.0 tells us that the given tank 
parameters are not consistent with a stabilizing tank. An 
example of greater than -2.0 is shown on figure 15 for the 
low tank on Ship no. 1. 

g. Tank Placement . There is another aspect of vertical 

tank placement in addition to comments already given about how 
it affects and GM. A special case of tank operation is 

whenever the duct is closed for some reason while the tank is 
full, for example for emergency maintenance. On high tanks 
especially there is a large GM loss due to the KG rise associated 
with the large tank fluid weight high in the ship. Figure 16A 
gives for an example the high tank on Ship No. 2 to show the 
undesirable effects of this operational condition. 

h. Ship Size . It seems obvious from the results to state 
that a lighter displacement ship such as Ship No. 1 apparently 

is not as well suited a candidate for an antirolling tank instal- 
lation as a larger ship. Whether this is due to the center of 
buoyancy being less submerged or to the inertia of the ship 
allowing for a higher natural frequency remains to be answered. 

7 . General Conclusions 

Although none of the design examples is complete in itself - 
primarily because a range of speeds and sea states was not looked 
at for each ship - it was possible to draw a number of general 
conclusions from them. 
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First of all, the equations of motion seem to be valid. 
Uncoupling sway by setting Xb = 0 and l/v* st 2 = 0 gives the clas- 
sical fourth order differential equation for a coupled ship-tank 
system. The assumed excitations are the primary departure from 
the latter, along with the sway coordinate, and hence introduce 
the greatest possible source of error because the reactive force 
and moment terms were also taken from classical theory. A lot 
of numerical accuracy hinges upon the two critical assumptions 
of uncoupling of roll and sway and application of the excita- 
tions at the center of (static) buoyancy. Model testing of some 
sort is needed to prove or disprove these conjectures. 

The general design method seems to be flexible and effi- 
cient within the range of interest. This was due to the veri- 

1/3 

fication that <)> versus y^. and K^., a 'drawable' three dimen- 
sional function, is a relatively flat function. This was due 
to its being an integral function. 

The solution of the fifth order differential equations of 
motion by the design method predicts tank parameter values which 
are quite similar and in some cases identical to the classical 
fourth order approach. The roll responses of the two approaches 
is not significantly different. However, a more refined solu- 
tion is possible by the methods of this paper. For example, the 
high tank of Ship No. 2 with y^ = 1.2 had a lower weight than 
that for the classically designed tank of y t = 0.97 due to the 
smaller cross duct associated with a higher tank frequency. 
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VII . Conclusions and Recommendations 
1. Conclusions 

The design method proposed in this paper seems to offer an 
attractive, efficient and flexible method for the design of pas- 
sive antirolling tanks. The analysis for active antirolling tanks 
led instead to a method for a feasibility study for the use of 
active tanks. The active tank design method predicts perform- 
ance of an economical active tank and will select optimum tank 
parameters for it. However, the design of the pump and the pump 
controller was left as a separate step. The active tank analysis 
also demonstrated why active tanks are such unpopular creatures. 
The small difference in tank angle which the pump is able to 
achieve and the small amount of active damping which it adds to 
the equations of motion are not nearly enough to overcome the 
chief cause of tank efficiency suffering in high sea states — 
the saturation of the tank which places an upper limit on the 
tank righting moment. 

The design method for passive antirolling tanks is attrac- 
tive because it uses a good representation of the real ocean 
wave process in predicting roll angle. Since there is much 
available information on ocean conditions for specific ship- 
ping routes or operational areas, this information may be sub- 
stituted for the assumed Pierson-Moskowitz function. The design 
method is efficient because the equations of motion lend them- 
selves to a complex algebraic form of solution, thereby resulting 
in a high numerical efficiency for the computer program. Because 
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of this a large range of solutions may be investigated at little 
cost. Finally, the design method is flexible in that tank de- 
sign may be approached from a wide variety of given requirements. 
For example, meeting some criteria for operational use, meeting 
space or weight constraints, or any combination of these may be 
imposed as design requirements. Also, the equations for the 
hydrodynamic and tank coefficients may be changed as necessary 
in the main computer program without altering the subroutine 
which calculates roll response. 

The prediction of roll angles seems reasonably well in agree- 
ment with all of the references, both in magnitude and in dis- 
tribution over heading angles. The short-crested sea predictions 
should be more accurate than those for long-crested or regular 
waves because linearity is better preserved. Only small compo- 
nents of roll for each heading angle are added together by super- 
position, including the worst heading at resonance peaks. 

A general conclusion concerning design parameters is that 
a reasonable passive tank design for preliminary design purposes 
may be estimated for a given X , or X = 0.20, u , from Appendix B, 

L. TH. L. 

from the Chadwick-Klotter equation on page 53, K^_ = -0.80, 
and the tank dimensions from Appendix B. 

2 . Recommendations 

It was decided to first investigate and validate the equa- 
tions of motion by taking the simplest case of a passive tank 
and linearized coefficients. As it became apparent that the re- 
sults of this approach were promising, a design method for passive 
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tanks was devised. A thorough checking of the results by model 
or full scale testing seems to be in order as the next step. 

The numerical integration technique seems to offer the best 
hope for modelling a first order, or higher order control effector 
in conjunction with the fifth order ship and tank equations. 
However, as previously mentioned this technique is inefficient 
in that there are so many integration steps to be performed in 
order to do an optimization, and several orders of magnitude 
more integration steps to be able to look at short crested seas 
over a range of sea states and ship speeds. In short, the 
numerical analysis technique has a certain merit for analysis 
but is a slumsy synthesis tool. Now that the basic equations 
for the reactive forces and moments have been validated, there 
may be a way to overcome this difficulty. This is to combine 
the spectral density function and the spreading function before- 
hand to make again an irregular, short-crested sea by taking 
all the harmonic frequencies and superposing them. A five or 
ten minute real time segment of excitations could be made and 
stored as library data. This data could be manipulated for 
speed, heading angle and sea state to form a simulation of the 
excitations. This could be used to drive the equations of motion 
with the effector included. Data for events such as significant 
roll angle could be obtained from the statistics of the results. 

It would be interesting to see how well the results of such a 
simulation compare with the complex algebraic solution. There 
is room for question as to whether the simulation would be as 
accurate as the complex algebraic solution because of the reduced 
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form some of the linear coefficients would assume. For example, 
the sway damping term must lose its frequency dependence in 
this type of approach because the excitation would be a series 
of absolute quantities derived from a frequency spectrum and 
would lack the frequency information necessary for calculation 
purposes . 

It was noticed when comparing the results of a passive tank 
design with the design solution contained in reference 13 that 
there was very close agreement in overall results in short- 
crested seas except for the case of the region of heading angles 
near 0° and near 180°; that is, head seas and stern seas. The 
author believes that the discrepancy most likely could be due 
to the simplistic approach taken towards obtaining a heading 
angle correction in Section III . 7. This author assumed a constant 
distribution of wave excitations over an effective shiplength. 

It is recommended that other distribution functions more closely 
resembling that which is actually found aboard the ship be tested 
to see whether a closer agreement could be realized. The impli- 
cations of using a second or third order parabolic distribution 
would be worth examining. Before doing this, however, it would 
be wise to carry through on the first recommendation given in 
this section and check on which of the predictions of roll angle 
are closer to model or full scale experimental tests, those of 
this author or those of reference 13. 

A final recommendation is that some of the steps for passive 
tank design could be better automated. For example, the inte- 
grals S' and S" could be accurately handled for any situation 
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by the proper numerical program. Also, the burden of solving 
the tank sizing equations on Figure 14 which need an iterative 
convergent solution should be removed from the tank designer. 
This would be no simple task of programming, however, because 
of the great number of combinations of constraints such as 
or w s ^. being limited by tank dimensional constraints. 
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Appendix A Anti-Rolling Tank Geometry Factors 

It was decided to construct a parametric family of Anti- 
Rolling tanks using the simple tank geometry of Figure A-l. In 
practice, a much more complicated geometry could conceiveably 
exist due to the need to use the ship's bulkheads, sides, and 
decks as part of the tank structure. Also, depending upon the 
tank damping required, it may be necessary to fair in sharp 
corners or other flow obstructions. It is expected, however, 
that the assumed geometry will be close enough to any required 
geometry in any specific installation so that equivalent dimen- 
sions and equivalent integrals (S' and S") can be quickly cal- 
culated and overall dimensions easily estimated. 

The basic geometry is shown below. 

I 




Figure A-l. Anti-Rolling Tank 
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The tank form factors are S' and S" . 

S' = q ^"(A °/a) dl 

S " = 0 5"( d / R ) dl 

Performing the required integrations by parts over each 
of the regions of constant sectional area, A q , , and A q , and 
dropping the nearly negligible higher order terms gives the 
desired results. 

S' = 2(L q + (R - b t )A Q /A 1 ) = 2R(L q /R + A q / A^ ) 

S" = 2R(2L /R + z /R) 
o t 

The two tank natural frequencies are u> t and to . 

•« t =V 2 g/ s ' 

« st =V2g/S" 

W t " V^7 r* * + Ao/Al) 

w st £ '/57T : yV(2Lo /R + 2 Vr) 

The above steps have reduced the problem of selecting a 
tank size to selecting a mix of the tank parameters R, (A /A^) , 
(L q /R) , and (z^/R). It may be seen that the other tank dimen- 
sions are known once these are established. The range of 
parameters for which this analysis is valid is approximately: 

0.2 < L q /R < 1.0 

- 1 < z t /R < + 2 

R < B/2 

Notice that for z^<-2L q the value of will be imaginary. 

Later work with the equations show that this indicates a 90° 
phase shift of the motion induced in the tank by the ship roll 



=Lo/R 
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Figure A-2 Tank Natural Frequencies 
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motion, but does not imply any dramatic change in the equations 
of motion. 

Figure A-2 summarizes the relationship between the tank 
frequencies and the tank parameters. For most combinations of 
these parameters, and in the range of validity previously men- 
tioned, Figure A-2 is sufficiently accurate for preliminary 
tank design purposes. 

A convenient way to determine tank size is by fluid weight. 

Tank weight (fluid) = 2pA Q R(L o /R + A^/A q ) 

An effectiveness paraneterof tank stabilizing potential is 
the ratio of tank righting moment to wave upsetting moment per 
unit of tank or waveslope angle. This ratio, A^., is also equal 
to the percentage reduction in static GM that the tank induces 
on the ship. 

At = 2pgA Q R 2 /K<}> = V 2 A q R 2 /35AGM (A in long tons, V = p/pH 2 0) 

Another measure of tank stabilizing potential is the maxi- 
mum waveslope that the tank could statically stabilize. Since 
the maximum tank angle possible is about L q /R (radians), the 
static tank capacity max is given by: 






max 



= AtL /R 
o 



rl> = 2) A L R/35AGM 
max o o 
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Appendix B. Coefficient Estimates 

1. Hydrodynamic Coefficients 

Figures B-l and B-2 display the hydrodynamic coefficients 
for ellipsoidal bodies of revolution, as theoretically pre- 
dicted by potential flow theory. The strategy in using these 
figures is to assume a body with the same dimensions (L,B,T,) 
as the ship. In cases where the mass of the body of revolution 
is quite a bit different than that of the actual ship, correction 
can be made by taking the average value of the coefficient for 
a body of the same dimensions (L,B,T) and for a body of the 
same mass as the ship with similar geometry (L/B, B/T, L/T) . 

A. Kp and Yv 

These are both constants. For order of magnitude 

estimates, Yv = -0.9m, and K. = -0.25 I . 

' P x 

B. K and Y 

_P X 

Kp, given as B44 on the figure, is assumed to be a 
constant mean value for reasons previously discussed. ^22 ° n 
the figure corresponds to Y ^ (io,U) by the notation of this 
paper. It is assumed to be well enough represented as a se- 
quence of linear segments. Figure B-3 graphically depicts this 
linear segmented approximation. The coordinates (F^, D^) are 
used as computer program inputs. The function is assumed 
constant after point 3; although this is certainly grossly 
inaccurate, the roll response operator for the normal range of 
ships is so small in this region that any assumption about the 
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behavior of the actual function will have essentially no effect 
on the roll angle solution. 

The sway damping component Y 2 (U ) is estimated by the 
methods of reference (3) to be: 

-irpT 2 U £ Yv2 (U) £ - (np/2) T 2 U 
The lesser absolute value is associated with a bare hull 
only, and the factor of two is added for a ship with a large 
amount of rudder surface area and other stabilizing underwater 
sail areas such as skegs, keels, etc. A particular ship may 
be somewhere between. For a normal amount of rudder area it 
will probably be close to the larger absolute value. 

C. K£ 

K<f> may be found directly from the formula K<J> = -mg GM. 
2. Ship Parameters 

A. tjy the center of buoyancy 

The center of buoyancy is accurately estimated by 
using Morrish's Formula for KB. 

t b = T - KB 

KB = 1/3 (5/2 T - 35A/Awp) 

Awp = Cwp • L • B 

Cwp = 0.237 + 0.833 Cp reference 14. 

B. GM and G. 

The vertical center of gravity and GM, if not known. 



may be estimated together by the relations: 
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KG = KB + BM - GM 
BM = Cit L (B) 3 /35A 

Cit = -0.013 + 0.10 Cp reference 14 

GM = 0.08 B (see reference 3 for a list of repre- 
sentative GMs ) 

Notice that if one or the other of G or GM are known, the 
unknown quantity is very accurately determined. Otherwise, 
both quantities suffer from being subject to 'guesstimates'. 

C. Ix 

The transverse moment of inertia of the mass of the 

~ 2 

ship may be closely estimated by: k = 0 . 38 B ; Ix = m (k ) 

X IK 
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Appendix C. Computer Program 
1. General Information 

A program which provides a numerical solution to the 
equations of motion was written in standard Fortran language. 
It is compatible with any Fortran capable system as long as 
the input device number KI and the output device number KO are 
correctly specified in the program. A plotting subroutine 
PICTR was used in the program; if this library subroutine is 
unavailable, replace the CALL PICTR cards with dummy CONTINUE 
cards. 

There are two basic programs. Main 1 and Main 2, 
which are separable decks and need not be run concurrently. 

A single subroutine RESP calculates the roll angle per unit 
waveslope, the tank angle per unit waveslope, and the phase 
angle between the tank and roll. A library subroutine, PICTR, 
is called to give a pictorial display, but because the results 
are printed out anyway, it is unessential to the program. 

Main 1 section of the program computes optimum tank 
damping and tuning for a given set of conditions. 

Main 2 section of the program computes roll response 
versus heading angle in short crested seas for a given ship, 
tank, and sea state. 

Because there are many untried variations to the 
program it is advised that anyone using it should proceed with 
care and check results at each step of the computations. 
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2. Computer Program Symbols 

The following symbol list cross references computer 
variables with the symbols used elsewhere in this text, where 
indicated. 

MAIN PROGRAM 



IX lx 

LB Ab 

KFI K<j> 

KPD Kp 

KY K 

y 

KS K 

s 



LY Ay 



KP 

H3 



K 

H 



P 

1/3 



LE Le 

C Spreading Function 

Coefficients 

XLAB Axis Labels 

TITLE Title Tables 

XMT jit 

XXMT Dummy y t 

F13 4> 1//3 

RAO Response Ampl. 
Operator 

0 a) 



S or SK iji 2/6u> 

AB Dummy holder 

E Tank Angle p 

TANG n 1//3 

KO Output device number 

KI Input device number 

BG BG 

GM GM 

CB Cb 

XL L 

B B 

T T 

TB t b 

XYVD | Yv/m | 

XLT At 

YV2 Y 

V 

U U 



XKT1 K t initial value 

DKTl increment 

XMST2 yst 2 
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MAIN PROGRAM (continued) 



DMT 




yt increment 




G2 


/ 2 


W1 




to initial value 


WS 


w s 


DW 




to increment, 


A to 


WE 


toe 


XKPD 




1 Kp/Ix | 




W 


too 


BETA 




3 




XI 


X 


D1 - 


D3 


Sway Damping 
values 


Y , 
vl 


EP 


Heading angle cor- 
rection, 1/2 X 


FI - 


F3 


to values correspon- 
ding to Di 


CTB 


-2tbto^ // g 
e ^ 


YV 




Yv 




ETA or E 


n 


YVD 




Y. 

v 




YV1 


Y 

V 


Gl 




U 




N 


Number of increments 
of A to 








SUBROUTINE 


RESP 




LT 




Xt 




W 


too 


LO 




Xo 




WS 


tos 


LY 




Xy 




Gl 


n 


LB 




Xb 




G2 


1 2 


M 




y 




TOCBL 


Phase angle 


M2 




2 

y 




Aij 


Matrix elements. 






o 






real 


MT2 




yt 
















Bij 


Matrix elements. 


MST2 




yst 






imaginary 


KY 




K 




Ci 


Intermediate matrix 






y 






steps 


KS 




K 

s 




RX 


<P/V 


KT 




K t 




RT 


T)/4> 



107 



3. Computer Program Input 

The following input data are required by the program. 
Ship Dimensions 



L 


Ship Length 




B 


Ship Beam 




T 


Ship Draft 




Cb 


Block Coefficient (Immersed volume/LBT) 




BG 


Distance from center of gravity to center 
ancy, + in z direction 


of buoy- 


GM 


Metacentric Height 




fc b 


Draft to the center of buoyancy 




Le 


Effective shiplength (between 60% beam sections) 




Tank Parameters 




yt 


Tank tuning (wt/ios) 




Ayt 


Tank tuning increment for search routine 




Kt 


Tank damping 




AK t 


Tank damping increment for search routine 




4-2 

yst 


Ship-tank decoupling frequency; primarily 
of vertical placement 


a function 


At 


Ship-tank coupling parameter -- % loss in static 
GM of ship due to tank free surface effect 



K 

P 

I | 

|Kj/Ix| 
Dl - D3 



Equat i on Coefficients 
Roll damping 

Added mass coefficient in sway. 
Added inertia coefficient in roll. 
Sway damping ordinates for Y v ^(to,U) 
Sway damping corresponding to 



FI - F3 
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Excitation Par a meters 

wl Initial value (lower cut-off frequency, about 0.2) 

Aw Increment of w (notice, 30 increments are used un- 

less N is changed in the program; the 30 increments 
plus w^ will determine the upper cut-off freq.) 

1/3 

H Significant waveheight of seas. 

8 Heading angle of ship with respect to the seas 

(180° for head seas) 

U Ship speed in the x-direction 

Program Control Variables 

There are three control variables which must be specified 
on the input deck. They are used to switch various functions 
in the program, and to terminate the computations. 

(+1, no PICTR) 

KSTOP ±1 Read Type 1 data cards (-1, PICTR at yt = ytl, 

Kt = Ktl for RAO for MAIN 1) 
(-1, PICTR of short crested 
seas for MAIN 2) 





0 


Stop 


execution of the program! ! 


JSTOP 


+1, 0 


Main 


1 Program executed 




-1 


Main 


2 Program executed 


LSTOP 


+ 1 


Read 


Type 2 data cards 




0 


Stop 


execution of program! ! 




-1 


Read Type 1 data cards; switch back to Main 
Program 



Examples of sample data decks will resolve any confusion 
caused by the above statements. The intent of having the 
switching variables was so as to be able to run Main 1, Main 2, 
or both together, and at the same time use some of the 
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subscripted computer variables (matrix elements) as dummy 
holders for intermediate calculation steps. Unless the data 
decks examples are followed, there are probably several non- 
workable combinations of the switching variables. 

Data Card Input Deck Format 

The general format for all the variables is floating point 
decimal. the exceptions are JSTOP, KSTOP, and LSTOP, which 
are fixed point decimal (12). 

There are three different types of data cards. To run 
Main 1, types 0, 1, & 3 are used. To run Main 2, types 0, 1, 

2, & 3 are needed. The data values on types 1 and 2 need not 
be identical; the reason for having type 2 data cards was to 
permit a number of variations of data without retyping a whole 
type 1 deck. 
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Card # 



Cplumn 



1 - 10 11 - 20 21 - 30 31 - 40 41 - 50 51 - 60 61 - 70 71 - 80 

(column 1-2 KSTOP) 

(column 3-4 JSTOP) 



1A 

IB 

1C 

ID 

IE 



BG 



K 



yst 



D1 



GM 


Cb 


L 


B 


T 




H 1/3 


1 | 


At 


blank 


Le 


U 


K ti 


AK t 


ytl 


A yt 


col 


Aco 


|k./i 1 

P X 1 


3 


blank 


D2 


D3 


FI 


F2 


F3 


blank 


blank 



/<* X-axis label, colums 1-40 — Y-axis label, columns 41-80 — 



2A 

2B 

2C 



H 



1/3 



yt 



K. 



(column 41-42 LSTOP) 43-80 blank 



/&- 



-Title to be printed below short crested sea results, col 1-80 — P/ 
-X-axis lable colums 1-40 — Y-axis label, columns 41-80 £>/ 






-tolank card- 



Combinations of data cards for successful computer runs: 

1. Find optimum tank parameters (Main 1) 

(1A - IF), (1A - IF), (1A - IF), 3; set KSTOP = 01, JSTOP = 00 

2. Find response versus heading angle (Main 2) 

(1A - IF) , (2 A - 2C) , (2A - 2C) (1A - IF) , (2A - 2C) (2A - 2C) ,3 

Set KSTOP = 01 or -1, JSTOP = 01, LSTOP = 01 or -1. 
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7. Sample Output Data 



SIG WAVE HT = 10.0 HEADING ANGLE =90.0 DEG 

SHIP SPEED = 43.00 EFFECTIVE SHIP LENGTH = 260.0. 

SHIP NAT FREQUENCY = 0.5142 SHIP-TANK COUPLING = 0.2000 
GM= 3.360 SHIP-TANK DECOUPLING FREQ= 1.04 
LB = 0.2167 LY= 8.9732 Gl= 3.3155 G2= 11.5796 

KS = 0.0453 KY=-2. 33446 MAS= 231884. IX= 75130336. 



KT 


FI1/3 


-0.400 


3.544 


-1.400 


5.783 


-0.900 


4.870 


-0.400 


3.544 


0.100 


4.966 


-0.150 


4.401 


-0.400 


3.544 


-0.650 


4.089 


-0.525 


3.733 


-0.400 


3.544 


-0.275 


3.706 


MT= 1. 


01 Fll/3= 


TANK LINEAR DAMPING 


KT 


FIl/3 


-0.400 


3.474 


-1.400 


4.202 


-0.900 
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-0.400 
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0.100 
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4.300 


-0.400 


3.474 


-0.650 


3.374 


-0.900 
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-0.775 


3.472 


-0.650 


3.374 


MT= 1. 


51 FIl/3= 


TANK LINEAR DAMPING 


KT 


FIl/3 


-0.400 


7.451 


-1.400 


6.939 


-2.400 


6.. 430 


-3.400 


6.190 


-4.400 


6.067 


-5.400 


5.996 


-6.400 


5.952 


-7.400 


5.922 


MT= 0. 


51 f 11/3- 



TANK LINEAR DAMPING = 
MI (MIN) + 1.35 FIl/3= 



3.71 ETA1/ 3= 5.67 

= -0.27 

2 



3.37 ETAl/3= 6.33 

= 0.65 

3 



5.92 ETAl/3= 0.32 

-7.40 

-1.8098 
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LONG CRESTED SEAS; SPEED= 



MT= KT= LT= 

1.20 - 1.00 0.20 

WS= LB= KS= 

0.40 0.22 -0.06 



HEADING ANGLE SIG ROLL ANGLE 



34.00 SIG WAVE HT= 15.0) 



SIG TANK ANGLE (DEGREES) 
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0.00 
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180.0 
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0.00 



SHORT CRESTED SEAS 



MT= 

1 . 20 
WS= 
0.40 



KT= 

- 1.00 

LB= 

0.22 



LT= 

0.20 

KS= 

-0.06 



HEADING ANGLE SIG ROLL ANGLE SIG TANK ANGLE (DEGREES) 
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SHIP # 2 



HIGH TANK 



LT=0 . 20 
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Appendix D. Other Solution Methods 

1 . Analog Solution 

The equations of motion were put into their proper non- 
dimensional form for electrical analog simulation. However, 
it was considered necessary to further investigate the numer- 
ical solution in order to properly normalize the non-linear 
variables, and especially coefficient expressions. Using the 
results of this paper, an analog simulation is a trivial 
problem if the modeling of the coefficients can be 
accomplished . 

2 . Numerical Integration 

In order to perform this technique, the three coupled 
equations forming a fifth-order system must be algebraically 
broken down into five coupled first order differential equa- 
tions. This was done using the set of dimensionless equa- 
tions in preference to the dimensional set. The five state 

variables chosen to represent the fifth-order system were: 

* • 

(f> , <}>, v* *, g , and n • One of the excitation terms, the iX Q 

term, was dropped to make the algebra a little more tractable. 
A parametric study was done in order to justify this step. 

The below figure shows that for the test case investigated 
for a "normal" sized ship, the iX Q term has practically no 
bearing on the excitation for w q > 0.2, which was previously 
shown to be the lower boundary of the excitation frequency 



Phase Angle 
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range of interest. Therefore, the sway excitation terra was 
shortened to X i p and roll excitation to (l-A.A ) ip . This kept 

y ^ y 

the excitation term uni-dimensional . 




Yv + Yv 

•U) U) | 



Yv 






Figure Dl. Sway Excitations. 



The five coupled ordinary differential equations, in 
matrix form, were found to be, 
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The above matrix equation is separated into a ship- tank 
linear matrix, and excitation matrix, a controller matrix, 
and a ship-tank non-linear matrix. The coefficients in the 
matrices are all constants for a - given set of conditions. A 
separate equation for p* should be added for a particular 
contoller . 

The identities for the 29 new coefficients, z^-z^, 

<5^-6 0 , and a^-a^ are given below. 



a 



1 




2 




a 



3 a 



a 



1 



2 





32 - “[1 + 





^5 X t* 1+a 3^ 
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B 4 X t K t 



^5 

6 5 = [a 16^ " Afc] 



a, a _X. X, y, X, X 

6 - = Al + X X + X (-5-| — + -^-^)] + -^4 

6 6 o b y y y w+ . 2 y 0 y 2 p„, 2 y . 



yt ' 2 



yt '2 



6_ = [a,/ + X. ] 6 0 = 



" a i k Q 



7 1 16 



8 6 



o 



a. 

S = —T 



B 

p o 



Vl 6 l X t 6 l 

Z 1 = [ ° n — - + x, K + 2 ] z „ — [ 



6 



o 



b Y B yt a 2 



X b Y l B 2 X t { 2 , 

2 '■ B ° P A, 

yt 2 



XvYiP, A ^ 6 -5 

Z 3 = - K y + -M- ] 

° v yt a 2 



z = r X b Y l B 4 X t 6 4 , 

z 4 1 6 2 J 

° V yt « 2 



X b Y l B 5 X t B 5 

z 5 = [ -nr-^ + -H“ ] 
° V a 2 



r X b Y l B 6 , , X t 6 6 , 

z 6 = l ~ 6 + X y 2— ] 

° V a 2 



r X b Y l B 7 , X t 6 7 , 
Z 7 = [ — 6 + T“ ] 

° V a 2 



Z 8 = [ ‘ 



' X b K Q Y l X t 6 8 1 
2 J 
y yt a 2 



6 , 



Z 9 = [ 



Vi (1 /V + Vi 



6 , 



y yt a 2 
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The matrix equations may be drawn as a signal flow graph 
in block form. The signal flow graph is useful for inter- 
preting exactly what the basic relationships are between all 
the variables. This information cannot be extracted from the 
simpler form of the three coupled equations. 




Figure D2. Ship Anti-Rolling Tank System 
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Simulation of roll by the use of the above equations 
without the necessity of directly using the spectral density 
function technique could be accomplished by taking an excita- 
tion made of randomly phased harmonic waves , in their proper 
proportions, and integrating for a ship time duration of 
enough cycles, say a thousand, to get a meaningful sample for 

determining d> . In the interests of amount of effort 
rms 

required, this should be done subsequently to finding the 
approximate solution by the complex algebraic method adopted 
by the author. 
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